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Abstract. 

The main purpose of this paper is to develop a unified approach of multi-parameter Hardy space 
theory using the discrete Littlewood-Paley-Stein analysis in the setting of implicit multi-parameter 
structure. It is motivated by the goal to establish and develop the Hardy space theory for the flag 
singular integral operators studied by Muller-Ricci-Stein [MRS] and Nagel-Ricci-Stein [NRS]. This 
approach enables us to avoid the use of transference method of Coifman- Weiss [CW] as often used 
in the LP theory for p > 1 and establish the Hardy spaces H^p and its dual spaces associated with 
the flag singular integral operators for all < p < 1. We also prove the boundedness of flag singular 
integral operators on BMOp and and from to for all < p < 1 without using the 

deep atomic decomposition. As a result, it bypasses the use of Journe's type covering lemma in this 
implicit multi-parameter structure. The method used here provides alternate approaches of those 
developed by Chang-R. Fefferman [CFl-3], Chang [Ch], R. Fefferman [F], Journe [Jl-2], Pipher [P] in 
their important work in pure product setting. A Calderon-Zygmund decomposition and interpolation 
theorem are also proved for the implicit multi-parameter Hardy spaces. 
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1. Introduction and statement of results 

The multi-parameter structures play a significant role in Fourier analysis. On the one hand, 
the classical Calderon-Zygmund theory can be regarded as centering around the Hardy-Littlewood 
maximal operator and certain singular integrals which commute with the usual dilations on R^, 
given hj d ■ X — {dxi, dxn) for d > 0. On the other hand, if we consider the multi-parameter 
dilations on i?", given hy 5 ■ x = {Six^, 5nXn), where 5 = {5i, 5„) e = then these 

n-parameter dilations are naturally associated with the strong maximal function ([JMZ]), given 

by 

(1-1) M,(/)(a;) = sup^ / \f{y)\dy, 

R 

where the supremum is taken over the family of all rectangles with sides parallel to the axes. 

This multi-parameter pure product theory has been developed by many authors over the 
past thirty years or so. For Calderon-Zygmund theory in this setting, one considers operators 
of the form Tf = K * f, where K is homogeneous, that is, 5i...5nK{5 ■ x) = K{x), or, more 
generally, K{x) satisfies the certain differential inequalities and cancellation conditions such that 
6i...6nK{d ■ x) also satisfy the same bounds. This type of operators has been the subject of 
extensive investigations in the literature, see for instances the fundamental works of Gundy-Stein 
([OS]), R. Fefferman and Stein [FSl], R. Fefferman ([F]), Chang and R. Fefferman ([CFl], [CF2], 
[CF3]), Journe ([Jl], [J2]), Pipher [P], etc. 

It is well-known that there is a basic obstacle to the pure product Hardy space theory and pure 
product BMO space. Indeed, the role of cubes in the classical atomic decomposition of H'p{R^) 
was replaced by arbitrary open sets of finite measures in the product H'p{R"' x R'^). Suggested 
by a counterexample constructed by L. Carleson [Car], the very deep product BMO{R^ x R^) 
and Hardy space Hp[R^ x R^) theory was developed by Chang and R. Fefferman ([Ch],[CF3]). 
Because of the complicated nature of atoms in product space, a certain geometric lemma, namely 
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Journe's covering lemma([Jl], [J2] and [P]), played an important role in the study of the bound- 
edness of product singular integrals on H'p{R'^ x K^) and BMO{K^ x K^). 

While great progress has been made in the case of pure product structure for both and 
HP theory, multi-parameter analysis has only been developed in recent years for the theory 
when the underlying multi -parameter structure is not explicit, but implicit, such as the flag 
multi-parameter structure studied in [MRS] and [NRS] . The main goal of this paper is to develop 
a theory of Hardy space in this setting. One of the main ideas of our program is to develop 
a discrete version of Caldcron reproducing formula associated with the given multiparameter 
structure, and thus prove a Plancherel-Polya type inequality in this setting. This discrete scheme 
of Littlewood-Paley-Stein analysis is particularly useful in dealing with the Hardy spaces for 
<p < 1. 

We now recall two instances of implicit multiparameter structures which are of interest to us 
in this paper. We begin with reviewing one of these cases first. In the work of MuUer-Ricci- Stein 
[MRS], by considering an implicit multi-parameter structure on Heisenberg(-type) groups, the 
Marcinkiewicz multipliers on the Heisenberg groups yield a new class of flag singular integrals. 
To be more precise, let m{C,iT) be the Marcinkiewicz multiplier operator, where C is the sub- 
Laplacian, T is the central element of the Lie algebra on the Heisenberg group H"^, and m satisfies 
the Marcinkiewicz conditions. It was proved in [MRS] that the kernel of m{£,,iT) satisfies the 
standard one-parameter Calderon-Zygmund type estimates associated with automorphic dilations 
in the region where \t\ < \z\'^, and the multi-parameter product kernel in the region where \t\ > \z\'^ 
on the space C"^ x R. The proof of the L^, 1 < p < oo, boundedness of m(£, iT) given in [MRS] 
requires lifting the operator to a larger group, x R. This lifts K, the kernel of m(£, iT) on 
EI"', to a product kernel K on H" x R. The lifted kernel K is constructed so that it projects to 
Khj 

oo 
— OO 

taken in the sense of distributions. 

The operator T corresponding to product kernel K can be dealt with in terms of tensor 
products of operators, and one can obtain their L^, 1 < p < oo, boundedness by the known pure 
product theory. Finally, the L^, 1 < p < oo, boundedness of operator with kernel K follows from 
transference method of Coifman and Weiss ([CW]), using the projection tt : H"^ x i? — > M" by 

n{{z, t), u) = {z,t + u). 

Another example of implicit multi-parameter structure is the fiag singular integrals on R^ x R^ 
studied by Nagel-Ricci-Stein [NRS]. The simplest form of fiag singular integral kernel K{x,y) on 
R^ x R^ is defined through a projection of a product kernel K{x, y, z) defined on x R^ 

given by 

(1.2) K{x,y)= J K{x,y-z,z)dz. 

A more general definition of fiag singular kernel was introduced in [NRS], see more details of 
definitions and applications of fiag singular integrals there. We will also briefiy recall them later 
in the introduction. Note that convolution with a fiag singular kernel is a special case of product 
singular kernel. As a consequence, the L^, 1 < p < oo, boundedness of fiag singular integral follows 
directly from the product theory on R^ x i?"^. We note the regularity satisfied by fiag singular 
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kernels is better than that of the product singular kernels. More precisely, the singularity of the 
standard pure product kernel on R^^xK^, is sets {{x, 0)}U{(0, y)} while the singularity of K{x, y), 
the flag singular kernel on x R^ defined by (1.2), is a flag set given by {(0, 0)} C {(0, y)}. For 
example, Ki{x, y) = is a product kernel on R^ and K2{x, y) — ^^^j^^^^ is a flag kernel on i?^. 

The work of [NRS] suggests that a satisfactory Hardy space theory shoifld be developed and 
boundedness of flag singular integrals on such spaces should be established. Thus some natural 
questions arise. Prom now on, we will use the subscript "F" to express function spaces or functions 
associated with the multi-parameter flag structure without further explanation. 

Question 1: What is the analogous estimate when p = 11 Namely, do we have a satisfactory 
flag Hardy space Hp{R^ x R^) theory associated with the flag singular integral operators? More 
generally, can we develop the flag Hardy space H^{R^ x R^) theory for all < p < 1 such that 
the flag singular integral operators are bounded on such spaces? 

Question 2: Do we have a boundedness result on a certain type of BMOf{R^ x R^) space 
for flag singular integral operators considered in [NRS]? Namely, does an endpoint estimate of 
the result by Nagel-Ricci-Stein hold when p = oo? 

Question 3: What is the duality theory of so deflned flag Hardy space? More precisely, do 
we have an analogue of BMO and Lipchitz type function spaces which are dual spaces of the flag 
Hardy spaces. 

Question 4: Is there a Calderon-Zygmund decomposition in terms of functions in flag Hardy 
spaces Hp{R^ x R^)7 Furthermore, is there a satisfactory theory of interpolation on such spaces? 

Question 5: What is the difference and relationship between the Hardy space H'p{R^ x R^) 
in the pure product setting and H^{R^ x R^) in flag multiparameter setting? 

The original goal of our work is to address these questions. As in the theory for p > 1 
considered in [MRS], one is naturally tempted to establish the Hardy space theory under the 
implicit multi-parameter structure associated with the flag singular kernel by lifting method to 
the pure product setting together with the transference method in [CW]. However, this direct 
lifting method is not adaptable directly to the case of p < 1 because the transference method is 
not known to be valid when p <1. This suggests that a different approach in dealing with the 
Hardy Hp{R'^ x R^) space associated with this implicit multi-parameter structure is necessary. 
This motivated our work in this paper. In fact, we will develop a unifled approach to study multi- 
parameter Hardy space theory. Our approach will be carried out in the order of the following 
steps. 

(1) We first establish the theory of Littlcwood-Palcy-Stein square fimction gj? associated with 
the implicit multi-parameter structure and the estimates of gp {1 < P < oo). We then develop 
a discrete Calderon reproducing formula and a Plancherel-Polya type inequality in a test function 
space associated to this structure. As in the classical case of pure product setting, these 
estimates can be used to provide a new proof of Nagel-Ricci-Stein's LP{1 < p < oo) boundedenss 
of flag singular integral operators. 

(2) We next develop the theory of Hardy spaces Hp associated to the multi-parameter flag 
structures and the bovmdcdness of flag singular integrals on these spaces; We then establish the 
boundedness of flag singular integrals from Hp to L^. We refer to the reader the work of product 
multi-parameter Hardy space theory by Chang-R. Fefferman [CFl-3], R. Fefferman [Fl-3], Journe 
[Jl-2] and Pipher [P]. 
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(3) We then establish the duality theory of the flag Hardy space and introduce the dual 
space CMO^p^ in particular, the duality of H]^ and the space BMOp- We then establish the 
boundedness of flag singular integrals on BMOp- It is worthwhile to point out that in the classical 
one-parameter or pure product case, BMO{R^) or BMO{R^ x R^) is related to the Carleson 
measure. The space CMOp for all < p < 1, as the dual space of Hp introduced in this paper, 
is then defined by a generalized Carleson measure. 

(4) We further establish a Calderon-Zygmund decomposition lemma for any H^{R^ x R^) 
function (0 < p < oo) in terms of functions in HpiR"" x R"^) and H^iR"^ x i?"^) with < 
Pi < p < P2 < oo. Then an interpolation theorem is established between H^i^K^ x R^) and 
Hf{R'^ X i?"^) for any < p2 < Pi < oo (it is noted that i?^(i?" x i?"^) = Lp{R'^+'^) for 
1 <p < oo). 

In the present paper, we will use the above approach to study the Hardy space theory associ- 
ated with the implicit multi-parameter structures induced by the flag singular integrals. We now 
describe our approach and results in more details. 

We flrst introduce the continuous version of the Littlewood-Paley- Stein square function gp- 
Inspired by the idea of lifting method of proving the L'p{R^ x R^) boundedness given in [MRS], 
we will use a lifting method to construct a test function defined on R^ x i?"^, given by the 
non-standard convolution *2 on the second variable only: 

,(1) =t„W,(2),'^ ,,\ ^ / J,W(^ n,- 



(1.3) '4^{x,y)='4^^'U2'^^'\x,y)= j ^^'\x,y - z)i^^'\z)dz, 

where V^^^ e 5(i?"+"^), V^^^ e 5(i?"^), and satisfy 

^|V^)(2-^ei,2-^6)l' = l 

3 

for all (a, ^2) e i?" X i?"^\{(0,0)}, and 

5^|^(^)(2-S)|^ = 1 

k 

for all Tj G -R"^\{0}, and the moment conditions 



Rn+m R-r 



for all multi-indices a, /3, and 7. We remark here that it is this subtle convolution *2 which provides 
a rich theory for the implicit multi-parameter analysis. 

For / e LP,1 < p < 00, gpif), the Littlewood-Paley-Stein square function of /, is defined by 



where functions 
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il:f\x,y) = 2("+^)^V^'^(2^a^,2^y) and ^^(z) = 2'^''i;^^\2'' z). 

We remark here that the terminology "imphcit multi-parameter structure" is clear from the 
fact that the dilation i/^j is not induced from 'i/;{x,y) explicitly. 

By taking the Fourier transform, it is easy to see the following continuous version of the 
Calderon reproducing formula holds on L^(_R"^"'""^), 

(1-5) f{x, y) = ^Yl ^J''^ * ^J''^ * 

j k 

Note that if one considers the summation on the right hand side of (1.5) as an operator then, by 
the construction of function i/^, it is a flag singular integral and has the implicit multi-parameter 
structure as mentioned before. Using iteration and the vector-valued Littlewood-Paley-Stein 
estimate together with the Calderon reproducing formula on allows us to obtain the L^, 1 < 
p < oo, estimates of gp. 

Theorem 1.1:. Let 1 < p < oo. Then there exist constants Ci and C2 depending on p such that 
for 

Ci\\f\\p<\\gF{mp<C2\\f\\p. 



In order to state our results for flag singular integrals, we need to recall some definitions given 
in [NRS]. Following closely from [NRS], we begin with the definitions of a class of distributions 
on an Euclidean space . A k — normalized bump function on a space is a C'^— function 
supported on the unit ball with C^— norm bounded by 1. As pointed out in [NRS], the definitions 
given below are independent of the choices of k, and thus we will simply refer to "normalized 
bump function" without specifying k. 

For the sake of simplicity of presentations, we will restrict our considerations to the case 
R^ = EP'+'^ X R^. We will rephrase Definition 2.1.1 in [NRS] of product kernel in this case as 
follows: 

Definition 1.2:. A product kernel on R^'^^^xR'^ is a distribution K on J?"^+"^+"^ njhich coincides 
with a C°° function away from the coordinate subspaces (0, 0, z) and (x, y, 0), where (0, 0) e 

and {x,y) G K^^"^ , and satisfies 

(1) (Differential Inequalities) For any multi-indices a = (ai, • • • , an), P = {Pi, ■ ■ ■ , Pm) o,nd 

Im = (71, • • • ,7m) 

\d^dld2K{x,y,z)\ < Co.,p,-y{\x\ + \y\)-n-m-\a\-m . |^|— -I7I 

for all {x,y,z) e R^ x x with \x\ + ^ and \z\ ^ 0. 

(2) (Cancellation Condition) 

I / d^d^yK{x,y,z)(t>i{5z)dz\ < + \y\)-n-m-\a\-\0\ 

for all multi-indices a,P and every normalized bump function 0i on R^ and every S > 0; 

I [ d2K{x,y,z)(l)2{Sx,6y)dxdy\<Cy\z\-'^-^^^ 

J 
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for every multi-index 7 and every normalized bump function 02 on and every S > 0; and 

K{x, y, z)(f)3{Six, Siy, S2z)dxdydz\ < C 



Ir 



' ^'n-\-m-\-r, 

for every normalized bump function 03 on and every Si > and ^2 > 0. 

Definition 1.3:. A flag kernel on x i?™ is a distribution on R^+'^ which coincides with a 
C°° function away from the coordinate subspace {(0, y)} C R^^'^ , where e R^ and y e R^ and 
satisfies 

(1) (Differential Inequalities) For any multi-indices a = (ai, ■ ■ ■ , an), P = (Pi, ■ ■ ■ , Pm) 

\d^d^K{x,y)\ < l«l -(1x1 + 1^1 

for all {x,y) e R"^ x R"" with \x\ ^ 0. 

(2) (Cancellation Condition) 



R" 



for every multi-index a and every normalized bump function 0i on R^ and every S > 0; 

I / d^K{x,y)MSx)dx\ < C^|y|---I^l 
Jr^ 

for every multi-index P and every normalized bump function 02 on i?" and every S > 0; and 

K{x,y)(f)3{Six,S2y)dxdy\ < C 



for every normalized bump function 03 on R^+"^ and every Si > and ^2 > 0. 

By a result in [MRS] , we may assume first that a flag kernel K lies in L^(i2"+™). Thus, there 
exists a product kernel K'^ on R^+'^ x R^ such that 



K{x,y) = / K'^{x,y — z, z)dz. 

J R^ 



Conversely, if a product kernel K'^ lies in L^{R^~^^ x i?"^), then K(xj y) defined as above is a flag 
kernel on R^ x R"^. As pointed out in [MRS], we may always assume that K[x, y), a flag kernel, 
is integrable on i?" x i?"^ by using a smooth truncation argument. 

As a consequence of Theorem 1.1, we give a new proof of the L^, 1 < p < 00, boundedness 
of flag singular integrals due to Nagel, Ricci and Stein in [NRS]. More precisely, let T{ f){x,y) = 
K * f(x, y) be a flag singular integral on i?" x R^. Then is a projection of a product kernel 
K« on X i?^. 
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Theorem 1.4:. Suppose that T is a flag singular integral defi,ned on x with the flag 
kernel K{x,y) = J K\x,y — z,z)dz, where the product kernel satisfies the conditions of 

Definition 1.2 above. Then T is hounded on for 1 < p < oo. Moreover, there exists a constant 
C depending on p such that for f & L^,! < p < oo, 

\\Tif)\\p<C\\f\\,. 

In order to use the Littlewood-Paley-Stein square function (/i? to define the Hardy space, one 
needs to extend the Littlewood-Paley-Stein square function to be defined on a suitable distribution 
space. For this purpose, we first introduce the product test function space on R^+'^ x R^. 

Definition 1.5:. A Schwartz test function f{x,y,z) defined on R^ x R^ x R^ is said to be a 
product test function on R^+'^ x R^ if 

(1.6) J f{x,y,z)x''y^dxdy = J f{x,y,z)z''dz = 

for all multi-indices a,(3,'-f of nonnegative integers. 

If f is a product test function on R^+'^ x R^ we denote f e Soo{R^^^ x R^) and the norm 
of f is defined by the norm of Schwartz test function. 

We now define the test function space Sp on R^ x R'^ associated with the fiag structure. 

Definition 1.6:. A function /(x, y) defined on R^ x R^ is said to be a test function in Sf{R^ x 
i?™') if there exists a function /* G S^oiR^^^ x R^) such that 

(1.7) f{x,y)= J f\x,y-z,z)dz. 

If f & Sf{R^ X R^), then the norm of f is defined by 

\\f\\sF{R"xR^) = infill/*^ II 5oo(i?"+-xi?-) : for all representations of f in (1.7)}. 
We denote by {Sp)' the dual space of Sp. 

We would like to point out that the implicit multi-parameter structure is involved in Sp- Since 
the functions ipj^k constructed above belong to Sp{R'^ x R^), so the Littlewood-Paley-Stein square 
function gp can be defined for all distributions in (Sp)'. Formally, we can define the flag Hardy 
space as follows. 

Definition 1.7:. LetO<p<l. i?^(i?" x i?"^) = {/ G (Sp)' : gp{f) G LP(i?" x i?^)} . 

/// G Hp(R^ X i?™-), the norm of f is defined by 

(1-8) ll/lk^ = \\9Fif)\\p. 

A natural question arises whether this definition is independent of the choice of functions if^j^k- 
Moreover, to study the i/^-boundedness of flag singular integrals and establish the duality result 
of Hp, this formal deflnition is not sufficiently good. We need to discretize the norm of Hp. In 
order to obtain such a discrete Hp norm we will prove the Plancherel-Polya-type inequalities. 
The main tool to provide such inequalities is the Calderon reproducing formula (1.5). To be more 
specific, we will prove that the formula (1.5) still holds on test function space Sp{R^ x R"^) and 
its dual space (Sp)' (see Theorem 3.6 below). Furthermore, using an approximation procedure 
and the almost orthogonality argument, we prove the following discrete Calderon reproducing 
formula. 
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Theorem 1.8:. Suppose that are the same as in (1-4)- Then 

(1-9) fix,y) = J^|/||J|V'j,fc(a;,y,a;/,?/j)V'j,fe*/(a;7,yj) 

j k J I 

where '^'j^ki^iViXiiyj) G Sf{R^ x R"^),I C R^,J C i?"^, are dyadic cubes with side-length 
£{I) = 2~^~^ and £{J) = 2~^~^ + 2~^~^ for a fixed large integer x/, yj are any fixed points 
in I, J, respectively, and the series in (1.9) converges in the norm of Sf{R^ x R^) and in the 
dual space [Sp)' ■ 

The discrete Calderon reproducing formula (1.9) provides the following Plancherel-Polya-type 
inequalities. We use the notation Ak. B io denote that two quantities A and B are comparable 
independent of other substantial quantities involved in the context. 

Theorem 1.9:. Suppose e 5(i?"+^), 0^^^ e 5(i?^) and 

ip{x,y)= J 'ip^^\x,y- z)'i/}^'^\z)dz, 



cf>{x,y)= J ct>^'\x,y-z)iP^^\z)dz, 



and il^jk, (f)jk satisfy the conditions in (1-4)- Then for f G (Sf)' and < p < oo, 



j k J I "e^.^e-^ 



(1-10) ^ llCEEE \<p,,k*f{u,v)Wj{x)xj{y)}'^\\, 

jkJI''^''''^' 

where ifjj^kix^y) and (j)j^k{x,y) are defined as in (1-4), I C R^,J C i?"^, are dyadic cubes with 
side-length £{I) = 2~^~^ and £{J) = 2~^~^ + 2~^~^ for a fixed large integer N,xi and xj are 
indicator functions of I and J, respectively. 

The Plancherel-Polya-type inequalities in Theorem 1.9 give the discrete Littlewood-Paley- 
Stein square function 



(1.11) g'^FifXx, y) = <! E E E E IV'i.fe * fi^i^ yj)\\i{^)xj{y) 

j k J I 

where /, J, xj, and yj are the same as in Theorem 1.8 and Theorem 1.9. 

From this it is easy to see that the Hardy space Hp in (1.8) is well defined and the Hp norm 
of / is equivalent to the norm of gp. By use of the Plancherel-Polya-type inequalities, we will 
prove the boundedness of flag singular integrals on H^. 
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Theorem 1.10:. Suppose that T is a flag singular integral with the kernel K{x, y) satisfying the 
same conditions as in Theorem I.4. Then T is bounded on Hp, for < p < 1. Namely, for all 
< p < 1 there exists a constant Cp such that 

m/)iiH^<Cpii/ii/f- 

To obtain the Hp boundedness of flag singular integrals, we prove the following general 

result: 



Theorem 1.11. Let 0<p<l. If T is a linear operator which is bounded on L^(i?"+'^) and 
HPp{W X i?"^), then T can be extended to a bounded operator from x K^) to LP(i?"+^). 

From the proof, we can see that this general result holds in a very broad setting, which 
includes the classical one-parameter and product Hardy spaces and the Hardy spaces on spaces 
of homogeneous type. 

In particular, for flag singular integral we can deduce from this general result the following 



Corollary 1.12:. Let T be a flag singular integral as in Theorem 1.4- Then T is bounded from 
H^p{W X R^) to LP(i?"+"") forQ<p<l. 

To study the duality of H^, we introduce the space CMOp. 

Definition 1.13:. Let ipj^k be the same as in (1-4)- We say that f e CMOp if f e {Sp)' and it 
has the finite norm \\f\\cMO'' defined by 



(1.12) 




for all open sets il in x with finite measures, and I C R^, J C R^, are dyadic cubes 
with side-length £{I) = and £{J) — + respectively. 



Note that the Carleson measure condition is used and the implicit multi-parameter structure 
is involved in CMO^p space. When p = 1, as usual, we denote by BMOp the space CMO\^. To 
see the space CMOp is well defined, one needs to show the definition of CMOp is independent 
of the choice of the functions V'i.fc- This can be proved, again as in the Hardy space Hp, by the 
following Plancherel-Polya-type inequality. 



Theorem 1.14:. Suppose ip,(f) satisfy the same conditions as in Theorem 1.9. Then for f e 



sup _ _ 



1 



(1.13) 
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where I C R^, J C i?"^, are dyadic cubes with side-length = 2~^~^ and i{J) = 2~^~^ + 
for a fixed large integer N respectively, and VL are all open sets in x R"^ with finite 
measures. 

To show that space CMOp is the dual space of Hp, we also need to introduce the sequence 
spaces. 

Definition 1.15:. Let be the collection of all sequences s — {sjxj} such that 



< oo, 



\s\\sp = <i^^\sixjf\i\ Vl ^xi{x)xj{y) 

j,k I, J 

LP 

where the sum runs over all dyadic cubes I C -R", J C R"^ with side-length £{I) — 2~^~^ and 
i{J) = 2~^~^ + 2~^~^ for a fixed large integer N, and x/, cind xj indicator functions of I 
and J respectively. 

Let cF be the collection of all sequences s = {sjxj} such that 
(1.14) ||s||^p=sup< — tzyX] X] k/xJp > < oo, 

" (M" j,k I,J:IXJCQ J 

where Q are all open sets in i?" x R^ with finite measures and the sum runs over all dyadic cubes 
/ C i?", J C i?"^, with side-length = 2'^'^ and 1{J) = 2'^'^ + 2-^-^ for a fixed large 
integer N . 

We would like to point out again that certain dyadic rectangles used in and reflect 
the implicit multi-parameter structure. Moreover, the Carleson measure condition is used in the 
definition of c^. Next, we obtain the following duality theorem. 

Theorem 1.16:. Let < p < 1. Then we have (s^)* = c^. More precisely, the map which maps 
s = {sixj} to < s,t >= X] sixjtixJ defines a continuous linear functional on with operator 

IxJ 

norm ||^||(sp)* ~ \\t\\cp, and moreover, every £ e (s^)* is of this form for some t G c^. 

When p = 1, this theorem in the one-parameter setting on R'^ was proved in [FJ]. The 
proof given in [FJ] depends on estimates of certain distribution functions, which seems to be 
difficult to apply to the multi-parameter case. For all < p < 1 we give a simple and more 
constructive proof of Theorem 1.16, which uses the stopping time argument for sequence spaces. 
Theorem 1.16 together with the discrete Calderon reproducing formula and the Plancherel-Polya- 
type inequalities yields the duality of Hp. 

Theorem 1.17:. LetO <p<l. Then (H^)* = CMO^p. More precisely, if g e CMO^p, the map 
£g given by £g{f) —< f,g >, defined initially for f e Sf, extends to a continuous linear functional 
on Hp with \\£g\\ ~ ||^||cmo^- Conversely, for every £ e (Hp)* there exists some g e CMOp so 
that £ = £g. In particular, {Hp)* = BMOp- 

As a consequence of the duality of Hp and the i^^-boundedness of flag singular integrals, we 
obtain the i?MOi?-boundedness of flag singular integrals. Furthermore, we will see that L°° C 
BMOp and, hence, the BMOp boundedness of flag singular integrals is also obtained. 

These provide the endpoint results of those in [MRS] and [NRS]. These can be summarized as 
follows: 
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Theorem 1.18:. Suppose that T is a flag singular integral as in Theorem 1.4- Then T is bounded 
on BMOf- Moreover, there exists a constant C such that 

mf)\\BMOF<C\\f\\BMOF. 

Next we prove the C alder on- Zygmund decomposition and interpolation theorems on the flag 
Hardy spaces. We note that x i?^) = for 1< p < oo. 

Theorem 1.19. (C alder on- Zygmund decomposition for flag Hardy spaces) Let < p2 < 1,^2 < 
p < Pi < oo and let a > be given and f G Hp{R'^ x R'^). Then we may write f = g + b 
where g G Hp{R'^ x i?"^) with p < pi < oo and b e Hp{R"' x i?"") with < p2 < p such that 
WdW^rlpi < C'cK^^"^! I/I l^p and ||fc||^P2 < C'ck^^^^I |/| |^p ? where C is an absolute constant. 

rip, p Hp Hp, 



Theorem 1.20. (Interpolation theorem on flag Hardy spaces) Let < p2 < pi < oo and T be 
a linear operator which is bounded from to L^^ and bounded from Hp to L^^ , then T is 
bounded from Hp to LP for all p2 < p < pi. Similarly, if T is bounded on Hp and Hp , then T 
is bounded on Hp for all p2 < p < pi. 

We point out that the Calderon-Zygmund decomposition in pure product domains was estab- 
lished for all L^ functions (1 < p < 2) into H^ and L^ functions by using atomic decomposition 
on iy-*^ space (see for more precise statement in Section 6). 

We end the introduction of this paper with the following remarks. First of all, our approach 
in this paper will enable us to revisit the pure product multi-parameter theory using the corre- 
sponding discrete Littlewood-Paley- Stein theory. This will provide alternative proofs of some of 
the known results of Chang-R. Fefferman, R. Fefferman, Journe, Pipher and establish some new 
results in the pure product setting. We will clarify all these in the future. Second, as we can 
see from the definition of flag kernels, the regularity satisfled by flag singular kernels is better 
than that of the product singular kernels. It is thus natural to conjecture that the Hardy space 
associated with flag singular integrals should be larger than the classical pure product Hardy 
space. This is indeed the case. In fact, if we deflne the flag kernel on i?" x i?"^ by 

K{x,y)^ J K{x - z,z,y)dz, 

where K{x,z,y) is a pure product kernel on R^ x i?"""""^, and let Hp be the flag Hardy space 
associated with this structure, thus we have shown in a forthcoming paper that HP{R^ x i?"*) = 

H^iR"" X R"^) n H^iR"^ X R"^). Results in [MRS] and [NRS] together with those in this paper 
demonstrate that the implicit multi-parameter structure, the geometric property of sets of singu- 
larities and regularities of singular kernels and multipliers are closely related. Third, the authors 
have carried out in [HL] the discrete Littlewood-Paley-Stein analysis and Hardy space theory in 
the multi-parameter structure induced by the Zygmund dilation and proved the endpoint esti- 
mates such as boundedness of singular integral operators considered by Ricci-Stein [RS] on H^ 
{0 < p < I) and BMOz, the Hardy and BMO spaces associated with the Zygmund dilation, e.g., 
on R\ given by 5{x, y, z) = {5ix, 52y, Si52z), di,52 > 0, where the L^ {1 < p < oo) boundedness 
has been established (see [RS] and [FP]). 

This paper is organized as follows. In Section 2, we establish the LP estimates for the multi- 
parameter Littlewood-Paley-Stein g— function for 1 < p < oo and prove Theorems 1.1 and 1.4. 
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In Section 3 we first introduce the test function spaces associated with the multi-parameter flag 
structure and show the Calderon reproducing formula in (1.5) still holds on test function space 
SpiR"' X K^) and its dual space {Sp)', and then prove the almost orthogonality estimates and 
establish the discrete Calderon reproducing formula on the test function spaces, i.e.. Theorem 1.8. 
Some crucial strong maximal function estimates are given (e.g. Lemma 3.7) and together with 
the discrete Calderon reproducing formula we derive the Plancherel-Polya-type inequalities, i.e.. 
Theorem 1.9. Section 4 deals with numerous properties of Hardy space Hp and a general result 
of bounding the norm of the function by its Hp norm (see Theorem 4.3), and then prove the 
Hp boundedness of flag singular integrals for all < p < 1, i.e.. Theorem 1.10. The boundedness 
from Hp to for all < p < 1 for the flag singular integral operators, i.e.. Theorem 1.11, is thus 
a consequence of Theorem 1.10 and Theorem 4.3. The duality of the Hardy space Hp is then 
established in Section 5. The boundedness of flag singular integral operators on BMOp space is 
also proved in Section 5. Thus the proofs of Theorem 1.14, 1.16, 1.17 and 1.18 will be all given in 
the Section 5. In Section 6, we prove a Calderon- Zygmund decomposition in flag multi-parameter 
setting and then derive an interpolation theorem. 

Acknowledgement. The authors wish to express their sincere thanks to Professor E. M. 
Stein for his encouragement over the past ten years to carry out the program of developing the 
Hardy space theory in the implicit multi-parameter structure and his suggestions during the 
course of this work. We also like to thank Professor J. Pipher for her interest in this work and 
her encouragement to us. 

2. estimates for Littlewood-Paley-Stein square function: Proofs of Theorems 1.1 
and 1.4 

The main purpose of this section is to show that the {p > 1) norm of / is equivalent to 
the norm of grif): and thus use this to provide a new proof of the boundedness of flag 
singular integral operators given in [MRS] . Our proof here is quite different from those in [MRS] 
in the sense that we do not need to apply the lifting procedure used in [MRS] directly. We flrst 
prove the estimate of the Littlewood-Paley-Stein square function gp- 

Proof of Theorem 1.1: The proof is similar to that in the pure product case given in [FS] 
and follows from iteration and standard vector-valued Littlcwood-Paley-Stein inequalities. To see 
this, deflne F : R^+^ H = P hy F{x,y) = {i/jj^^ * f{x, y)} with the norm 

\\F\\H = {J2\i^^'^*fix,y)\Y^. 
j 

When X is flxed, set 

9iF){x,y) = {J2U^k^ *2F(a;,.)(y)||l,}i 

k 

It is then easy to see that g{F){x, y) = gF{f){x, y). If x is flxed, by the vector-valued Littlewood- 
Paley-Stein inequality, 

J g{FY{x,y)dy<C j 

Jim Jim 

However, ||F||^ = IV']^'' * /(a;, y)P}^, so integrating with respect to x together with the 

i 
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standard Littlewood-Paley-Stein inequality yields 

j j gF{fnx,y)dydx<C J J {J^li^f * fi^,y)\V dydx < C\\frp, 

R" R" R<^ ^ 

which shows that ||^f(/)||p < 

The proof of the estimate < C\\gFU)\\v routine and it follows from the Calderon 

reproducing formula (1.5) on L2(i?"+^), for all / e n , e n L^' and i + ^ = 1, and 
the inequality ||5rF(/)||p < which was just proved. This completes the proof of Theorem 

1.1. Q.E.D. 

Remark 2.1: Let V^^^ e S{W+'^) be supported in the unit ball in and V^^^ e S{W^) be 

supported in the unit ball of and satisfy 

IV'«(^6,^6)l'^ = i 







for all (Ci,6) e i?" X i?"^\{(0,0)}, and 



Jo ^ 



for all?7 G R^\{<d]. We define V'*(a;,y, z) = ilj^^\x,y- z)ip^^\z). Set ijr{x,y) = ^-""^^^^^1 : f ) 
and 'ilji^\z) = s~"*V(f ) and 

^t,s{x,y)^ / ^lji^\x,y - z)^ljP{z)dz. 



I Jim 

Repeating the same proof as that of Theorem 1.1, we can get for 1 < p < oo 

oo oo 

(2.1) Uj l\i^t,s*f{x,y)\'j^}'^\\,<C\\f\\„ 



and 

oo oo 

(2.2) ll/l|p^ll{/ l\^t,s*^t,s*f{x,y)fj^}-^\\,. 





The LP boundedness of flag singular integrals is then a consequence of Theorem 1.1 and 
Remark 2.1. We give a detailed proof of this below. 

Proof of Theorem 1.4: We may first assume that K is integrable function and shall prove the 
L^, 1 < ]3 < oo, boundedness of T is independent of the norm of K. The conclusion for general 
K then follows by the argument used in [MRS]. For aU / e L*', 1 < p < oo, by (2.2) 

oo oo 

(2.3) \\T{mp<C\\{J J \i;,,,*i;,,,*K*f{x,y)\'j^}'^\\,. 
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Now we claim the following estimate: for / e L^, 1 < p < oo, 

(2.4) IV't.s * K * /(x, y)\ < CMs{f){x, y), 

where C is a constant which is independent of the norm of K and Mg{f) is the maximal 
function of / defined in the first section. 

Assuming (2.4) for the moment, we obtain from (2.3) 

oo oo 

(2.5) \\Tf\\, < C\\{ 1 1 (M,(V'M * DMfj^j'^Wp < C\\f\\„ 



where the last inequality follows from the Fefferman-Stein vector-valued maximal function and 
Remark 2.1. 

We now prove the claim (2.4). Note that 'ipt,s * K{x,y) = J ipf^g * K'^{x,y — z,z)dz, where 
1^1 g(x, y, z) is given in Remark 2.1 and K{x, u) = f K\x, y—z, z)dz, where K'^{x, y, z) is a product 
kernel satisfying the conditions of definition 2.1.1 in [NRS] (or Definition 1.2 in our paper). The 
estimate in (2.4) will follow by integrating with respect to z variable from the following estimate: 

(2-6) \iPlg*K\x,y,z)\<C- , , ^, „ ^ 7 , 

where the constant C is independent of the norm of K. The estimate (2.6) follows from that 
in the pure product setting R^+^ x given by R. Fefferman and Stein [FS]. Q.E.D. 

3. Test function spaces, almost orthogonality estimates and discrete Calderon repro- 
ducing formula: Proofs of Theorems 1.8 and 1.9 

In this section, we develop the discrete Calderon reproducing formula and the Plancherel- 
Polya-type inequalities on test function spaces. These are crucial tools in establishing the theory 
of Hardy spaces associated with the fiag type multi-parameter dilation structure. The key ideas 
to provide the discrete Calderon reproducing formula and the Plancherel-Polya-type inequalities 
are the continuous version of the Calderon reproducing formula on test function spaces and the 
almost orthogonality estimates. 

To be more precise, we say that a function a(a;, y) for (a;, y) G x R^ belongs to the class 
S°°{R^ X R^) if a{x,y) is smooth and satisfies the difierential inequalities 



(3.1) md^a{x, y)\ < ^iv,a,/3(l + \x - y\ 



-N 



y 

and the cancellation conditions 



(3.2) J a{x,y)x°'dx = J a{x,y)y^dy = 

for all positive integers N and multi-indices a,/? of nonnegative integers. 

The following almost orthogonality estimate is the simplest one and its proof can be adapted 
to the more complicated orthogonal estimates in subsequent steps. 
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Lemma 3.1. If ip and (p are in the class iS°°(i?"' x R^), then for any given positive integers L 
and M, there exists a constant C — C{L,M) depending only on L,M and the constants Ajv,a,/3 
in (3.1) such that for all t, s > 

(3.3) I j M^, ^)M^, y)dz\ < C{- A ^ '^"^ 



s t' (tVs + |a;-?/|)("+-^)' 

where iptix, z) = t~"'V(f > f ) ^'^^ 0s (-2, y) = •s~'^0(f )? and t As = min{t, s),t\/ s = max{t, s). 
Proof: We only consider the case M = L = 1 and t > s. Then 

/ il^t{x,z)4)s{z,y)dz 

= [^t{x,z) -il)t{x,y)](j)s{z,y)dz 

= [ +[ =1 + 11 

We use the smoothness condition for ipt and size condition for 0^ to estimate term /. To estimate 
term //, we use the size condition for both i/j-t and 0^. 

For the case M > 1 and L > 1, we only need to use the Taylor expansion of iptix,-) at y and 
use the moment condition of ipf We shall not give the details. Q.E.D. 

Similarly, if 'ip^{x, y, z, u, v, w) for (a;, y, 2;), (tt, v, w) E R'^ x R^ x R^ is a smooth function and 
satisfies the differential inequalities 

\d^^d^^d2^d:^d^^d2j^ijHx, y, z, u, V, w)\ 

(3.4) < ^iV,M,ai,a2,/3i,/32,7i,72(l + \x - u\ + \y - v])'^ {1 + \z - w])'^ 

and the cancellation conditions 

/ ^'(.. y, u, V, dxdj = / ^«(.. „, u, V, 

(3.5) = /^«(., y, „, ^.iu'^v^dndv = / ^'(^, y, „, .^H. = 0, 

and for fixed xq e -R", yo G R^, (f)'^{x, y, z, xq, yo) e SooiR^'^'^ x R^) and satisfies 

\d^^d^^d2^(P^{x,y,z,xo,yo)\ 



(3.6) < BM,M,a^,0^,'y^,{'^ + ^ " ^o| + Iv - yo|)"'^(l + \z\)-^ 

for all positive integers N, M and multi- indices cii, 0:2, /^i, /92, 7i, 72 of nonnegative integers. Then 
we have the following almost orthogonality estimate: 
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Lemma 3.2. For any given positive integers Li,L2 and Ki,K2, there exists a constant C = 
C{Li, L2, Ki, K2) depending only on Li, L2, Ki, K2 and the constants in (3.4) and (3.6) such 
that for all positive t, s, t' , s' we have 



J V'm(^) '^)^t',s' yo)dudvdw\ 



^'n-\-m-\-m 

(3.7) 

- r ^s' s' {tyt' + \x-xo\ + \y-yo\Y''+'^+^^Hsys' + \z\)('^+^^y 

where i;l{x,y, z,u,v,w) = t-"s-™V"(f , f , f , f , f , f ) and 

4^,ix,y,z,xo,yo) = |, ^, ^, ^, ). 

The proofs of the almost orthogonahty estimate in (3.7) is similar to that in (3.3). We will 
only provide a brief proof. 

Proof of Lemma 3.2: We only consider the case Li — L2 — Ki — K2 — 1, t > t' and s < s'. 
Thus 



j V'mI^^ 2/5 ^5 '^^ '^5 ''^)(Pt',s' i'^^ ^0, yo)dudvdw\ 

/ 



i?" X X it'" 

^ • Bdudvdw 

i?" X i?'" X i?'" 



where 
and 



^ = i^lsi^i 'f^' ^) - i^isi^i y^ ^o, yo, w) 



B = (j)l^^,{u,v,w,xo,yo) - (pl,^^,{u,v,z,xo,yo). 
In the above, we have used the cancelation properties 

J {u, V, w, xq, yo)u'^v^dudv = 0, J ipt,s{^j y-> ^1 '^■> '^^ w)w^ = 

fin X Rm Jim 

for all multi- indices a, j3 and 7. 
Next, 

I j A ■ B dudvdw 

Jln X Rm X Rm 

A ■ B dudvdw 



u-xo\ + \'"-yo\<^it+\x-xo\ + \y-yo\), \w-z\<^{s'+\z\) 



+ I A ■ B dudvdw 

-xo\ + \v-yo\<^it+\x-xo\ + \y-yo\), \w-z\>^{s'+\z\) 



f 

J lu—xol + lv- 

f 

J\u 

j 

J \u — 3 



+ I A ■ B dudvdw 

\u-xo\ + \v-yo\>^it+\x-xo\ + [y-yo[), \w-z[<^{s'+\z\) 



+ I A ■ B dudvdw 

'\u-xo\ + \y-yo\>hit+\x-xo\ + \y-yo\), \w-z\>^{s' + \z\) 

--I + II + III + IV 
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For term /, we use the smoothness conditions on ipf^g and 0^, g,; For term //, we use the smooth- 
ness conditions on ipl^g and size conditions on 0^, g, ; For term we use the size conditions on 
ipl^g and smoothness conditions on 0^, g,; For term IV, use the size conditions on both ipl^^ and 
0^, g,. We shall not provide the details here. Q.E.D. 

The crucial feature, however, is that the almost orthogonality estimate still holds for functions 
in Sf{R^ X R^). To see this, we first derive the relationship of convolutions on x R^ and 
^n+m ^ ^j^^ usc this relationship frequently in this paper. 

Lemma 3.3. Let ip,(pe SpiR"" x R""), and ip^, 0« e 5oo(i?''+"' x i?"^) such that 



Then 



/ 'ilj\x,y - z,z)dz, (f){x,y) ^ (f)\x,y - z, z)dz. 

J J 



Lemma 3.3 can be proved very easily. Using this lemma and the almost orthogonality estimates 
on X i?"^, we can get the following 

Lemma 3.4. For any given positive integers Li,L2 and Ki,K2, there exists a constant C = 
C{Li, L2, Ki, K2) depending only on Li, L2, Ki, K2 such that iftyf < s V s', then 

\il^t,s*<l>t',s'{x,y)\ 

t\r ,s s\r (tytY' isy s')^^ 



't' t' ^s' s' (tVt'+ |a;|)("+-^i) (sVs' + 
and if t V t' > s V s' , then 

\i^t,s *<Pt',s'{x,y)\ 

t\r ,8 S\r (tVt')^l UytY^ 



t' t' ^s' s' Vt' + (tVt' + 



Proof of Lemma 3.4: We first remark that we will prove this lemma with K2, L2 
replaced by K^, -C'i, L'2. Thus, we are given any fixed K[,K'2i L[,L'2. 

Note that 

i^t,s*(i>t',8'{x,y)= J 'iljlg*(t)l,^g,ix,y- z,z)dz, 

where G 5oo(i?"+"' x i?™), and 

V'm * = j ijls{x-u,y -v,z-w)(l)l,^g,{u,v,w)dudvdw, 

Jim X Rm X Rn 
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Then by the estimate in (3.7), for any given positive integers Li,L2 and Ki, K2, there exists 
constant C = C{Li, L2, Ki, K2) depending only on Li, L2, i^i, K2 such that 

\il^t,s*(l>t',s'{x,y)\ 

(3.8) 



- H' ^s' J {tVt'+\x\ + \y-z\)('^+^+^^) {sVs'+\z\)(^+^^) ' 

Case 1: If t V < s V s' and |y| > s V s' , write 



J {tVt'+ \x\ + \y- ^|)(n+^+-f^i) (s V s' + |^|)(^+-^2) 



It is easy to see that 



+ / =/ + // 

z\<h\yUor \z\>2\y\ |j/|<|2|<2|s/| 



<c 



{tVf + \x\ + \y\){n+m+Ki) 

(tVt')^^ (sVs')"'+'^2 



(iVi' + |a;|)('^+^i) (s V s' + |?/|)'^+^2 

where we have taken Ki = K[ + K2 and used the fact that tV t' < s V s' and > s V s'. 
Next, we estimate 

i„i <. / ..d. 



(s V S' + |y|)(^+-^2) J {tVf +\x\ + \y - 2|)(n+m+i^i) 

(sVs')^^ (tVt')^' 



<c 



(S V S' + \y\)i^+K2) (tyt' + |x|)('^+^i) 
(sVs')^^ (tVt')^^ 



where we have used K2 > K2 and Ki = K[ + K2> K[. 
Case 2: If t V < s V s' and < s V s', then 

(tVt')^' (sVs')^' 



{tyt' + \x\ + \y- z\)^'^+'^+^^) (s V s' + |2|)("^+-^2) 
it'" 

~ (s V s')"" 7 (t V t' + |a;| + |y - ^|)("+'^+-'^i) 



<C 



{s V gQ^^ (t V tQ^i 

(s V s' + |y|)^+-^2 (t V f + |a;|)("+-^i) ' 
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Case 3: We now consider the case tV t' > sV s' and \y\ < tV t'. Then 



/ 



{tVt'+ \x\ + \y- (s V s' + |^|)(^+-^2) 



dz 



< 



{tyt' + |a;|)(^+'^i) {tyt' + \x\ + ||/|)("^+^2) ■ 

by noticing that Ki = K[ + 

Case 4: If we assume tV t' > s\/ s' and \y\ >t\/t', then we divide the integral into two parts 
/ and // as in the case t\/ t' < sV s'. Thus we have 



<C 



{tVf + \x\ + \y\){n+m+Ki) 



where we have taken Ki = K[ + K^- 
To estimate //, we have 



\ii\ <c 
<c 



Q.E.D. 



Roughly speaking, il^t,s * (l>t' ,s' (x , y) satisfies the one-parameter almost orthogonality when 
tV t' > sV s' and the product multi-parameter almost orthogonality when t\/ t' < s V s'. More 
precisely, we have the following 

Corollary 3.5. Given any positive integers Li,L2,Ki and K2, there exists a constant C = 
C(Li,L2,Ki,K2) > such that 

(i) If t > s we obtain the one-parameter almost orthogonality 

and, ift<s the product multi-parameter almost orthogonality is given by 

(3.9) lA..>M.,y)\ < C(j, A -)-.(j7 A -)^.^_^^^_^^. 

(a) Similarly, ift' > s', 

t t' s s' t'^^ t'^^ 

(3.10) \iht s * (f)t' s'(x,y)\ < C(- A A -)^'7 , , ^ , , , ^ , , 
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and iff < s', 

(3.11) 1^,,. * < C(i A A 7)-' (,^|';|U..,(,,^|^^.K,) 



Corollary 3.5 is actually what we will use frequently in the subsequent parts of the paper. The 
proof of Corollary 3.5 is a case by case study and can be checked with patience. We shall omit 
the details of the proof here. All these estimates will be used to prove the following continuous 
version of the Calderon reproducing formula on test function space Sf{R^ x K^) and its dual 
space {Sf)'. 

Theorem 3.6. Suppose that i/jj^k o,re the same as in (1-4)- Then 

(3.12) f{x, y) = ^Yl ^J''^ * ^J''^ * 

j k 

where the series converges in the norm of Sf o,nd in dual space (Sf)'- 

Proof: Suppose f e Sf and f{x, y) ^ J f\x, y - z, z)dz, where e 5oo(-R""^"* x K^). Then, 
by the classical Calderon reproducing formula as mentioned in the first section, for all f'^ e L'^, 

(3.13) fKx, y,z) = ^Yl Ak * ^ik * fK^, v. A. 

j k 

where z) = iljf\x, y)'^f\z). 

We claim that the above series in (3.13) converges in S^oiR^^^ x R^). This claim yields 

Wfi^^y)- ^j,k*'^j,k* f{x,y)\\sp 

-N<j<N -M<k<M 

=\\ j [fK^^y- z,z) - Y Y '^],k*'^],k* fK^^y- ^^^)]dz\\sF 

-N<j<N -M<k<M 

<\\f\x,y,z)- Yl E ^ife*V'ife*/«(^,y,^)ll5^ 

-N<j<N -M<k<M 

where the last term above goes to zero as N and M tend to infinity by the above claim. 
To show the claim, it suffices to prove that all the following three summations 

E E4.*V'i.*/^E E E Y^l*4k*f^ 

|j|>Ar|fc|<M \j\<N\k\>M |i|>Ar|A;|>M 

tend to zero in Soo{R^~^"^ x R^) as N and M tend to infinity. Since all proofs are similar for 
each of the summations, we only prove the assertion for the first summation which we denote by 
/at m- Note that 

/1,m(^'?/'^) = E E / V'*,fe*^*,fe(a;-'">y-'f^>^-'^)/n«>'f^>'^)c^«c?^c?w 
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where tjj'^*ip'^{x~u,y — v,z — w) satisfies the conditions (3.4) and (3.5), and /'^(tt, w, to) satisfies the 
conditions (3.6) with xq = yo = 0. The almost orthogonahty estimate in (3.7) with t = 2~\s = 
2~^, t' = s' — 1 and xq — yo — 0, imphes 



/ 



'ipjj^*'ipjj^{x — u,y — v,z — w)f\u, V, w)dudvdw\ 

i?" X it'" X 



V 1)^1 



{2-^ V 1)^2 



{2-j V 1 + |a;| + \y\){n+m+K^) (2"^ V 1 + \z\f^+^^) 



This, by taking Li> Ki and L2 > ii'2, gives us 

lim sup (1 + + |y|)-+-+^^(l + k|)"^+^1/i m(^.2/, ^)I = 0- 



Since d'^d^d2{f\[ j^){x,y, z) = {d'^d^d^ f^)N,M{x,y, z) and applying the above estimate to 
d^d^djp which also satisfies the conditions in (3.6) with = j/o = 0, we obtain 

lim sup (1 + \x\ + |y|)-+-+^ni + \z\r+'''mdld2f^)N,M{^^y^^)\ = 0, 

N,M^oo xeR"^ ,yeR'^ ,zeR^ 

which shows the claim. 

The convergence in dual space follows from the duality argument. The proof of Theorem 3.6 
is complete. Q.E.D. 

Using Theorem 3.6, we prove the discrete Calderon reproducing formula. 

Proof of Theorem 1.8: We first discretize (3.12) as follows. For / e 5^, by (3.12) and using 
an idea similar to that of decomposition of the identity operator due to Coifman, we can rewrite 



f{x^ y) = Yl^ j J -u,y-w) (ipj^k * f) (w, w)dudw 

j,k I, J J J 



(3.14) 



=5:e 

j,k I, J 



J I 



{'(pj,k * /) {xi, yj) + n{f){x, y). 



We shall show that TZ is bounded on Sp with the small norm when / and J are dyadic cubes 
in i?" and with side length 2~^~^ and 2~^~^ + 2~^~^ for a large given integer and a;/, yj 
are any fixed points in /, J, respectively. 

To do this, assuming f{x,y) = J f^{x,y - z,z)dz, where /« e 5oo(i?"+'" x i?^). When 
k < j, we write 

ii^3,k * /) {u, w) - {ipj^k * f) {xi, yj) 
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where the last integral above is over i?" x x x K^. 
When k > j, we write 

{^j,k * f) {u, w) - (ijjj^k * f) {xi, yj) 

f^^ {u,w - v,v)dv - J (^i;^-,^* f^^ {xi,yj - v,v)dv 

We now set 



= ^ ^ / / ^j,k{x -u,y-w) [{ipj^k * f) {u, w) - {ipj^k * f) {xi, yj)] dudw 
j,k i,J J I 




TZ^ {x, y — z,z,u,v,w )y{u ,v ,w )du dv dw dz 
n\p){x,y-z,z)dz, 

where TZ^{x^ y, z, u', v', w') is the kernel of 7^* and 

7?.* {x, y — z, z,u' , v', w') 

= XU^/ / / i^j^\x -u,y - z -w)ij''^\z) 
k<j I, J J J 

x — u' ,v — v')ij/^\w — V — w') — ipj^\xi — u' ,v — v')ip^\yj — w' — v) dudwdv 

+ E E E E /77 ^^^'^ - ^' ^ - ^ - 

j k>j J I J I 

X j^V'j^'' {u — u' ,w — V — v') — il^^p (xi — u' ,yj — V — v') ip^^^ {v — w')dudwdv. 

Using the change of variables from ^ to ^ + v — in the term oik < j, and ^ to 2; — v in the term 
of /c > J, we can rewrite 

7?.* {x, y, z, u', v', w') 
= EE / / / ^f\x-u,y-v)%^f\z + v-w) 

k<j I,J j j 

X — u' ,v — v')'tpj^\w — w' — v) — ipj^\xi — u',v — v')ipj^\yj — w' — v) 

j k>j J I J I 

X |^'0j''"-*(it — u' ,w — V — v') — ipj^\xi — u' ,yj — V — i(^^\v — w')dudwdv 

= ^^jk+^Bjk- 

k<3 k>j 



dudwdv 
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We claim that 7^* is bounded in S^{R'^+"' x R"^). 
To see this, write 



J I 



J I 



X [ipj^^ {u — u' ,v — v') — {xi — u' ,v — v')]tp^^ {w — w' — v)dudwdv 
+ / / / t^f\x-u,y-v)%l)f\z + v-w) 



J I 



J I 



{xj ~ u' ,v — v')[ijj^^\w — w' — v) — i^^^\yj — w' — v)]dudwdv 
It is not difficult to check that ifj^^^ (z + v — w)'4>^^^ (w — w' — v)dw satisfies all the conditions 

(2) 

as V'fc (-2 — w') does with the comparable constants of Sf{R^ x R^) norm and that 

J J J i^j^\x~u,y — v)['iljj^\u — u',v — v') — 'iljj^\xi — u',v — v')]du^ 
^ J I 

satisfies all conditions as ifjj^^x — u'.y — v') but with the constants of SpiR^ x R^) norm replaced 
by C2^^ . This follows from the smoothness condition on il^^p (say the mean- value theorem) and 
the fact that u,xi and 1{I) = 2~^~K 

If we write 



J / 



J I 



dudv 



X / '0^^'* (2 + t> — w)iIj^^ {w — w' — v)dw 

J 

then the function A^^^{x, y, z, u', v\ w') satisfies all conditions as 

'i{jj^\x — u' ,y — v') ■ ip^^\z — w') 

does but with the Sf{R^ x R^) norm constant replaced by C2~^ . 
By the proof of Theorem 3.6, we conclude that 

J Af,^^ {x, y, z, u', v', w')f\u\ v', w')du'dv'dw' 

is a test function in Soo{R^^^ x R^) and its test function norm is bounded by C2~^ . 
Similarly, 

J + v- w)['il;l^\w - w' -v) - '4;f\yj - w' - v)]dw 
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satisfies all conditions as '4'^\z — w') but with the constant replaced by C2 ^ , which follows 
from the smoothness condition on ■^^^^ and the fact w,yj e J and 1{J) = 2~^~^ + 2~^~^ and 
k<j. 



We also note that 



j j j il^^p {x — u,y — v)ip^p {xi — u' ,v — v')dudv 



J I 



J I 



satisfies the same conditions as ^p^p {x — u',y — v') does with comparable Sf{R^ x EJ^) norm 
constant. Thus, we conclude that 

Af^ {x, y, z, u\ v', w')f^{u', v', w')du'dv'dw' 
is a test function in Sqiq{R"'~^^ X i?"^) and its test function norm is bounded by C2 ^ . Therefore, 

is a test function in Soo{R^~^^ x R^) and its test function norm is bounded by C2~^ . 
Similarly, we can conclude that 

Bjk{x, y, z, u', v', w')f\u', v', w')du'dv'dw' 

is also a test function in Soo{R^~^^ x R"^) and its test function norm is bounded by C2~^ . 
This shows that 7^«(/«)(a;, y, ^) e Soo{R''^'^ x R"") and 

which implies that 7^(/) e SpiR"" x R"") and 

(3.15) ||7^(/)||5.(i^"xi^-) < C2-^||/||5^(^„xfl-)- 

By (3.14) together with the boundedness of 71 on Sp with the norm at most C2~^ , if N is 
chosen large enough, then we obtain 



j k J I 



i=0 



^^7?.* / / i/jj^ki' ~ u, ■ — v)dudv 



J I 



Set 



i=0 



^^7?.* / / V'j,fe(' — u,- ~ v)dudv 



J I 



{x,y) = \I\\J\'i{;j,k{x,y,xi,yj). 



It remains to show il!j^kix,y,xi,yj) e Sp- This, however, follows easily from (3.15). 
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Q.E.D 

We next establish a relationship between test function in /* e <Soo (-R""'''^ x R'^) and test 
function f{x, y) = J f'^{x, y—z, z)dz under the actions oiTZ^^TZ and the implicit multi-parameter 

dilations. 



To be more precise, we define 



where 



Note that 



ni{{f%k)ix,y,z) 

R\x, y, z, u\ v\ w'){f)jk{u', v', w')du'dv'dw' 

2{n+m)j2mk J J J'R}(^X,y, z, u' , v' , w') f (2^ u' ,2^ v' , 2''w')du' dv' dw' 
n^{x, y, z, 2"%', 2-h', 2" v', w')du'dv'dw' 




and 

n\2-^x, 2-'y, 2-^z, 2-^u', 2-^v', 2" V) = 2("+™)^2'^'=7e»(a;, y, z, u', v' , w'). 
Thus we have 

n^{{P)j,k){x,y,z) = {n^{f%k{x,y,z). 

This implies 

n{fj,k){x,y) = f n\if^)j,k){x,y-z,z)dz 

= I {n^{f%k{x,y-z,z)dz={n{f))jk{x,y) 

It is worthwhile to point out that 

i^j,k{x, y, XI, yj) = j i^lkix, y- z,z, xi, yj)dz, 

where 

^,ki^,y,z,xi,yj) = {ip^)j,k{x,y,z,xi,yj), i)\x,y,z,xi,yj) e SooiBJ"^"^ x R^) 
and satisfies the condition in (3.6) with xq = xj, yo = yj. 
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Remark 3.1: If we begin with discretizing (3.12) by 

i J / ;{ 

and repeating the similar proof, then the discrete Calderon reproducing formula can also be given 
by the following form 

fix, y) = J2J2J2J2 l^W'^l'^jM^ -xi,y- yj)^j^f.{f){xi, yj), 
j k J I 

where \I\\J\'iljj^i^{f){xi, yj) = / / '4^j,k * ij^y{.f){'^i v)dudv. We leave the details of these proofs 

i=0 J / 

to the reader. 

Before we prove the Plancherel-Polya-type inequality, we first prove the following lemma. 

Lemma 3.7. Let /, J, J' be dyadic cubes in and respectively such that £{I) — 2~^~^ , 
e{J) = + 2-^-^, £{!') = 2-^'-^ and e{J') = 2"^'"^ + 2"'='-^. Thus for any u e I and 



2-i\J-f\ + \k-k'\)K2-ij'+k')K^j,njn 

S (2-r + lu - ...ir" (2-^' + 1« - ,,-1)"^'^ * 

<C . 2-l''-''l'f I M. (EEI'*^'.'''*/(^''.!'-'')Ix^'X/' ) I M 

k'<j' [ \ J' I' J } 

and 

2-i\j-fl + \''-k'\)K2-2j'Kipiljfl 

fe'>i'/',j'(2 ^ (2 ^ +\v-yj'\) 

where M is the Hardy-Littlewood maximal function on 7?"+™, Mg is the strong maximal function 
on R^ X R^ as defined in (1.1), and max |^^^, | < r. 

Proof: We set 

A, = {I' : = 2-^'-^, ^^-1^ < 1} 



and for £ > 1, z > 1 



2-3' 

Bo = W : £(J') = 2-^'-^ + 2-^'-^, < 1} 



Ae = {I' : £{!') = 2-^"-^, 2^-^ < < 2^}. 
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Bi = U' : £(J') = 2-^'-^ + 2-^'-^, T-^ < ^^4^ < 2n. 

2 



Then 

2-0"+fe')^|/'||j'| 



\4>j',k' *f{xi',yj')\ 



i2-3' + |« - xr\T^'' {2-k' + 1^ - yj,\r^'' 
< ^ 2-^(-+^)2-^(-+^)2-^("+-) * f{xr,yj 



e,i>0 I'£Ai,J'eBi 
£,i>0 \ \l'eAe,J'&Bi I 



<C{N) {mAJ2 \<l>j',k'*f{xi',yj'Wxi'XJ' {u,v) 



The last inequahty follows from the assumption that r > ^^j^ and r > which can be done 

by choosing K large enough. 

Similarly, we can prove the second inequality of the lemma. Q.E.D. 

We now are ready to give the 

Proof of Theorem 1.9: By Theorem 1.8, f e Sp can be represented by 

fix,y) = Y'Y'Y'Y\J'\\I'\^j',k'{x,y,xr,yj') {4>j'^k> * f) {xi',yj'). 
j' k' J' r 

We write 

{i^j,k * f) {u,v) 

= YYYY\I'\\J'\ (i^j,k*^j',k'i-,-,xi',yj'j^ iu,v){(f)j'^k' * f) (xr ,yj'). 
j' k' J' r 

By the almost orthogonality estimates in (3.10) and (3.11), for any given positive integer 
we have if j' > k' , 

\iljj,k*(^j',k'i-r)iu,v)\ 

2-j'K 2-^'^ 



(3.15) ^c2-^3-3'\K .2-\k-k'\^ 



(2-r + \u- xr\)'^+^ (2-^' + 1^; - yj'\)'^+^ 
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and ii j' < k', we have 

\ipj,k*^j',k'{-,-){u,v)\ 

O-j'K 0-j'K 

(3.16) <c2-b-^"l^.2-l'=-'='l^- 



(2-i' + |'u-a;//|)^+'^ + \v - yj>\)-^+^' 

Using Lemma 3.7, for any w, x// e /, f , yj/ e J, 
\'4^j,k* f{u,v)\ 



fc'<i' /' J' 



+ \u- xr\)''+^ {2-^' + \v- yj'\)'^+^ 

+^ E EE2"''"''"'2-l'=-'='l''|/'||J'|x 
k'>j' J' r 

2-j'K 2-^'^ 



\(t>j',k' *f{xr,yj' 



\4'j',k' * fixi',yj' 



k'<j' [ J' r ) 

+C 5] 2-l^-^-'l^2-l'=-'='l^|M('j]J]|0,.,fc,*/(x,^2/jOlx/'XJ'l I M 

k'>j' [ \ J' I' J } 

where M is the Hardy-Littlewood maximal function on i?""*""^, Mg is the strong maximal function 
on R"^ X K^, and max{^^^, ^j^ } < r < p. 

Applying the Holder's inequality and summing over j, k, I, J yields 



{EEEE \^j,k* f{u,v)Wixj 

[ j k J I "e/,^eJ 

[ j' k' I J' I' 

Since xj' and j/j/ are arbitrary points in /' and J', respectively, we have 



^EEEE \i^3,k* fiu,v)\'^xixj 

( j k J I «e-^.«eJ 



<^ { EE ^^(EE.J;-f,,, \<l>r,k' * /(«,^)|X/'X.')^ 

j' k' K J' I' 
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and hence, by the FefFerman-Stein vector-valued maximal function inequality [FS] with r < p, we 
get 



j k J I «e^.«eJ 



j' k' J' I' 



IP- 



This ends the proof of Theorem 1.9. Q.E.D. 

4. Discrete Littlewood-Paley-Stein square function, boundedness of flag singular 
integrals on Hardy spaces H^, from Hp to L^: Proofs of Theorems 1.10 and 1.11 

The main purpose of this section is to establish the Hardy space theory associated with the 
flag multi-parameter structure using the results we have proved in Section 3. As a consequence 
of Theorem 1.9, it is easy to see that the Hardy space Hp, is independent of the choice of the 
functions ip. Moreover, we have the following characterization of Hp using the discrete norm. 

Proposition 4.1. Let < p < 1. Then we have 

(4.1) ^ II I EEEEi^i-'^*/(^^'2/^)i'x/(^)xj(y) 

( j k J I 

where j,k,'4>iXiiXJi^iiyj ^'^^ same as in Theorem 1.9. 

Before we give the proof of the boundedness of flag singular integrals on Hp^ we show several 
properties of H^. 

Proposition 4.2. Sf{R^ x K^) is dense in Hp. 

Proof: Suppose / G H^, and set W = {(j, k, I, J) : \j\ < L,\k\ < M, I x J C S(0, r)}, where /, J 
are dyadic cubes in R^,R^ with side length 2~^~^ ,2~^~^ + 2~^~-^ , respectively, and -8(0, r) 
are balls in R^+'^ centered at the origin with radius r. It is easy to see that 

E l-^ll'^l xi, yj)^j,k * fixi, Vj) 
{j,k,i,j)ew 

is a test function in Sf{R'^ x R'^) for any fixed L,M,r. To show the proposition, it suffices to 
prove 

E \i\\J\i^j,ki^iyiXiiyj)^j,k* fixi,yj) 

ij,k,i,j)ew- 

tends to zero in the Hp norm as L,M,r tend to infinity. This follows from (4.1) and a similar 
proof as in the proof of Theorem 1.9. In fact, repeating the same proof as in theorem 1.9 yields 

II E \i\Wj,k{^^y^^i^yj)'^j,k* f{xi,yj)\\H^^ 

{j,k,i,J)ew- 
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<c\\{ Yl IV'i'fe * /(^^' yj)\^xixj}-' \\p, 

{j,k,I,J)&W- 

where the last term tends to zero as L, M, r tend to infinity whenever / e Hp. Q.E.D. 

As a consequence of Proposition 4.2, L'^{R'^~^^) is dense in Hp{R'^ x i?"^). Furthermore, we 
have 

Theorem 4.3. If f e L2(i?"+"^) n H^iR"^ x R""), < p < 1, then f e LP{R^+'^) and there is a 
constant Cp > which is independent of the norm of f such that 

(4.3) \\f\\p<C\\f\U^. 

To show theorem 4.3, we need a discrete Calderon reproducing formula on L^(i?"^"'""^). To be 
more precise, take (j)^^^ e CQ°(i?"^+'") with 

/ y)x'^y^dxdy = 0, for aUa, /3 satisfying < \a\ < Mq, < < Mq, 

J Jln+m 

where Mq is a large positive integer which will be determined later, and 

1^(2-^^1,2-^6)1' = 1, for aU(6,6) e i?"+^\{(0, 0)}, 

i 

and take 0(2) e C^{R^) with 

0^2) (2)2^^2 = for aU < I7I < Mq, 

and Efc W^K'^~^^2? = 1 for aU ^ e i?^\{0}. 

Furthermore, we may assume that (f)^^^ and 0^^) are radial functions and supported in the unit 
balls of R^+'^ and R^ respectively. Set again 



(f>jk{x,y)= / (j)f\x,y- z)(f)^^\z)dz. 

J 



By taking the Fourier transform, it is easy to see the following continuous version of Calderon 
reproducing formula on L^: for / e L2(i?"+"^), 

y) = YY ^jk * 0jfc * f{x, y). 

j k 

For our purpose, we need the discrete version of the above reproducing formula. 
Theorem 4.4. There exist functions (j)jk and an operator T^^ such that 

fix, y) WjM^ -xi^y- yj)4>i,k * {T^^if)) {xi, yj) 
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where functions (j)jk{x — xj, y — yj) satisfy the conditions in (3.6) with ai, 71, A^, M depending 
on Mq, xq = xj and yo = yj- Moreover, is bounded on L^(i?"'+'") and Hp{R^ x i?"^), and 

the series converges in L'^{R^~^"^). 

Remark 4.1: The difference between Theorem 4.4 and Theorem 1.8 are that our (f)jk in 
Theorem 4.4 has compact support. The price we pay here is that (pjk only satisfies the moment 
condition of finite order, unhke that in Theorem 1.8 where the moment condition of infinite order 
is satisfied. Moreover, the formula in Theorem 4.4 only holds on L^(i?""''™') while the formula in 
Theorem 1.18 holds in test function space Sp and its dual space {Sp)'- 

Proof of Theorem 4.4: Following the proof of Theorem 1.8, we have 
(4.4) fix,y) = J2J2J2J2^ / <Pj,k{x-'^^^y-^)dudv]{(pj,k* f)ixi,yj) + n{f){x,y). 



j k J I 



J I 



where /, J, j, k and TZ are the same as in Theorem 1.8. 
We need the following 

Lemma 4.5. Let < p < 1. Then the operator TZ is hounded on L^(i?"'+'^) fl H^p{R^ x W^) 
whenever Mq is chosen to he a large positive integer. Moreover, there exists a constant C > 
such that 

|W)||2<C2-^||/||2 

and 

\\^{f)\\Hl{R^xR^) <C2~ ||/||//|,(i?"xi?-)- 

Proof of Lemma 4.5: Following the proofs of Theorems 1.8 and 1.9 and using the discrete 
Calderon reproducing formula for / e L^(i?"'+™), we have 

\\9Fmf))\\p 



y 3 k J I 



E E l-^'ll^'ll (V'i.fc*^ ('0j',fe'(-,a://,-,2/j/) ■V'i'fc' * f{xi',yj')y^ ?XiXj 



j ^k,J,I j' ,k^ ,J\I^ 

where j, k, xi, XJi ^Ii yj ^i'^ the same as in Theorem 1.9. 
We claim: 



(i^j,k *T^{i^j',k'i-,xr,-,yj')^^ {x,y)\ 



0-{3M')K n-(kAk')K 

^^-N \K2-\k-k'\K _ / ^ ' ^ ^ ^~ 

(2-0-Aj') + \x-xj,\ + \y-z- yj,)n+m+K (2-(feAfe') + |^|)m+K 



R" 



where K < MQ.,max{:^^^, u^k ) ^ -^0 is chosen to be a lager integer later. 
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Assuming the claim for the moment, repeating a similar proof in Lemma 3.7 and then Theorem 
1.9, we obtain 



2 - 1 II 



j' k' J' I' 

j' k' J' I' 

It is clear that the above estimates still hold when p is replaced by 2. These imply the assertion 
of Lemma 4.5. 

We now prove the Claim. Again, by the proof of Theorem 1.8, 

^ [i^j',k' {■, xr, •, yj')) {x, y) = / TZ^{x, y- z,z, u', v', w')i)y^k' (•, a;//, •, yj')du'dv'dw'dz 

where 7?.^(x, y, z, u\ v\ w') is similar to TZ'^ as given in the proof of Theorem 1.8 but, as we pointed 
out in Remark 4.1, that the difference between TZ"^ here and Ti} given in the proof of Theorem 1.8 is 
the moment conditions. However, the almost orthogonality estimate still holds if we only require 
sufficiently high order of moment conditions. More precisely, if we replace the moment conditions 
in (3.5) "for aU ai, 71, q;2, /32, 72" by "for aU |q;i|, |7i|, \a2\, |/32|, I72I < Mq where Mq is a 
large integer, then the estimate in (3.7) still holds with Li, L2, i^i, K2 depending on Mq. Thus, 
the claim follows by applying the same proof as that of Theorem 1.8, and the proof of Lemma 
4.5 is complete. Q. E. D. 

We now return to the proof of Theorem 4.4. 

Denote {Tn)~^ = E»^i^*> where 

TNif) = ^^^^ij j ^3,k{x -u,y- v)dudvd] {^j^k * /) {xi, yj). 



j k J I 



J I 



Lemma 4.5 shows that if N is large enough, then both of Tn and (Tjv) ^ are bounded on 
L^(i?'^+'^) n H^{R^ X R^). Hence, we can get the following reproducing formula 

/(a;, y) = ^^^Y^ \I\\J\4>j,k{x -xi,y- yj)4>j,k * (^j^^/)) {xi, yj) 
j k J I 

where ^jkix - xi,y - yj) = -^^tjt J J (f)jk{x - xi - (u - xi),y - yj - {v - yj))dudv satisfies the 

J I 

estimate in (3.6) and the series converges in LP'^R^^'^). 

This completes the proof of Theorem 4.4. Q.E.D. 
As a consequence of Theorem 4.4, we obtain the following 

Corollary 4.6. /// e L2(i?'^+"^) n iy^(i?" x i?"^) andQ<p<l, then 

^ll{(E E E E \^3^ * ^^n'U)) i^i, yjrxiix)xj{y)}f^ \\p 

j k J I 
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where the constants are independent of the norm of f. 

To see the proof of Corollary 4.6, note that if / G L'^{R^~^"^), one can apply the Calderoh 
reproducing formulas in Theorem 1.8 and 4.4 and then repeat the same proof as in Theorem 1.9. 
We leave the details to the reader. We now start the 

Proof of Theorem 4.3: We define a square function by 

9{f){^,y) = {EEEE l'^^-.'^ * (T^'if)) {xi,yj)\\i{x)xj{y)V^ 
j k J I 

where (f)jk are the same as in Theorem 4.4. By CoroUary 4.6, for / e L'^{W+'^) n x W^) 

we have, 

<c\\f\\ 

To complete the proof of Theorem 4.3, let / e L^{R''+'^) n x i?"^). Set 

= {{x,y) eR^xR^: g{f){x,y) > T}. 

Denote 

Bi = {{j,k,i,j):\{ixj)nni\ > ^\ixJi\{ixj)nni+i\ < ^|/x j|}, 

where I,, J arc dyadic cubes in R^^R^ with side length 2~^~-^ ^2~^~^ + 2~^~^ , respectively. 
Since / e L^(i?"^+"^), by the discrete Calderon reproducing formula in Theorem 4.4, 

f{x,y) 

= E E E E ~xi,y- yj)\I\\J\(j)j,k * {T-\f)) {xi, yj) 

j k J I 

= E E \i\\J\^jA^-^i:y-yj)(p3,k* iTN^{f)){xi,yj), 

i {j,k,I,J)eBi 

where the series converges in norm, and hence it also converges almost everywhere. 
We claim 

II E l^WjM^ -xi,y- yj)cf>j,k * {xi^yjX < Cr^l^il 

{j,k,I,J)€Bi 

which together with the fact < p < 1 yields 

11/11^ < E II E l^ll<^l^i,'^(^ -xi,y- yMi,k * {T^\f)) {xi.yjWp 

i (j,k,LJ)eBi 

<CE2'^l^i| 

i 

<c\mw, < c\\f\f^,^. 

To show the claim, note that (p^^^ and i/'*-^-' are radial functions supported in unit balls. Hence, 
if {j, k, I, J) e Bi then (j)j^k{x — xj, y — yj) are supported in 



= {{x,y) : Ms{xQ,){x,y) > -^}. 
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Thus, by Holder's inequality, 

II Yl \J\mAx-xi,y-yj)<l>j,k*{T^\f)) {xi,yj)\\l 

{j,k,i,j)eBi 

<|n;|i-f|| \J\\I\kk{x-xi,y-yj)(l>j,k*{T^\f))ixj,yj)r^. 
{j,k,i,J)eBi 

By the duality argument, for all g E L'^ with ||^||2 < 1, 

l< \J\\I\^j,h{x-xi,y-yj)(j)j,k*{T^^{f)){xi,yj),g>\ 

{3,k,i,j)eB, 

= 1 Y \J\\I\^3,k*9{xi,yj)(j)j,k*{T^^{f)){xi,yj)\ 
ij,k,i,J)eBi 



2 / \ 2 



<c\ Y mmj,k*{TN\f))i^i,yj)\'] ■ Y \i\m,,k*9{xi,yj)\ 



2 



(j,k,i,j)eBi / \{3,k,i,J)eBi 



Since 



2 



< 



2 



Y \I\\J\\ct>j,k^g{xi,yj)\'' 
^{j,k,i,j)eBi 

Y \i\\J\ {m, (^^j,k*g) {x,y)xi{x)xj{y) 

Sj,k,I,J)&Bi 

1 

{j^Lr. Lr. (^^■''^ *^)' i^^y)d^dy)j < C\\9\\2 

thus the claim now follows from the fact that < C\Qi\ and the following estimate: 
C2^'\ni\> J g\f){x,y)dxdy 

{j,k,i,j)eBi 

>l Y \i\mj,k*{TN\f))i^i,yj)\^ 

(j,k,I,J)&Bi 

where the fact that |(/ x J) fl Oi\Oi+i| > i|/ x J| when {j,k,I,J) e Bi is used in the last 
inequality. This finishes the proof of Theorem 4.3. Q.E.D. 

As a consequence of Theorem 4.3, we have the following 

Corollary 4.7 x i?"^) is a subspace of Li(i?" x R'^). 
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Proof: Given / G Hp{R^^'^), by Proposition 4.2, there is a sequence {fn} such that fn G 
L2(i?"+"") n and converges to / in the norm of Hl{R^+'^). By Theorem 4.3, 

converges to g in L^(i?'^+"^) for some g G L^(i?"+"'). Therefore, / = in {Sp)' ■ Q.E.D. 

We now turn to the 

Proof of Theorem 1.10: We assume that K is the kernel of T. Applying the discrete Calderon 
reproducing formula in Theorem 4.4 implies that for / G L^(i?'^+'") n H^{R^ x i?*"), 

iliE E E E i'^^-'^ * ^ * /(^' v)\"xii^)xj{v)yn\p = 

j k J I 

II {E E E E I E E E E I i^U-.'^ * ^ * (• - ' • - ) y) X 

j k J I j' k' J' I' 

(l>j',k' * (TN^f)) ixi',yj')\'^xi{x)xAy)}^\\p, 

where the discrete Calderon reproducing formula in L^(i?"^+"^) is used. 

Note that (f)jk are dilations of bump functions, by estimates similar to the those in (2.5), one 
can easily check that 

10,- * K * 0,v,fc,(- -xr,-- yj'){x, y)\ < C2-l^-^'l^2-l'=-^'l^ 

J (2-0' ^i') + \x-xr\ + \y- z- yj'\Y+'^+^ ' + \z\)'^+^'^^' 

where K depends on Mq given in Theorem 4.4 and Mq is chosen to be large enough. Repeating 
a similar proof in Theorem 1.9 together with Corollary 4.6, we obtain 

\\Tf\\Hi < cii{EE{^^(EE I'^i'^'^' * i^N\f)) {^r.yj')\xj'xi'y}h^,y)}Hp 

j' k' J' r 

^ C\\{Y.Y.Y.Y.\^r,k' * [T-\f)) {xv.yj,)\\Ay)xi'{x)]H, < C\\f\U^, 
j' k' J' r 

where the last inequality follows from Corollary 4.6. 

Since L^(i?'^+'^) is dense in Hp{R^ x i?"^), T can extend to a bounded operator on H^[R^ x 
R"^). This ends the proof of Theorem 1.10. 

Proof of Theorem 1.11 We note that n is dense in i/^, so we only have to show this 
for / G H^p n L?. Thus Theorem 1.11 follows from Theorems 4.3 and 1.10 immediately. Q.E.D. 

5. Duahty of Hardy spaces and boundedness of flag singular integrals on BMOp 
space: Proofs of Theorems 1.14, 1.16 and 1.18 

This section deals with the duality theory of flag Hardy spaces Hp{R^ x R^) for all < p < 1. 
We first prove Theorem 1.14, the Plancherel-Polya inequalities for CMO^ space. 

Proof of Theorem 1.14: The idea of the proof of this theorem is, as in the proof of Theorem 1.9, 
again to use the discrete Calderon reproducing formula and the almost orthogonality estimate. 



DISCRETE LITTLEWOOD-PALEY-STEIN MULTI-PARAMETER ANALYSIS 



37 



For the reader's convenience we choose to present the proof of Theorem 1.14 in the case when 
n = m = 1. However, it will be clear from the proof that its extension to general n and m is 
straightforward. Moreover, to simplify notation, we denote fj^k — fn when R — I x J G and 
|/| = , \ J\ = 2~^~^ + 2~^~^ are dyadic intervals respectively. Here N is the same as in 

Theorem 1.8. We also denote by dist{I, I') the distance between intervals / and 

Sr= sup \i;R* f{u,v)\^, Tr= inf \(Pr * f{u,v)\^. 

u£l,veJ u£l,v£J 



With these notations, we can rewrite the discrete Calderon reproducing formula in (1.9) by 

fix,y)= Yl \A\J\4>R{^^y)(l>R* fi^i^yj)^ 



R=IxJ 



where the sum runs over all rectangles R = I x J. 

Let R' = r X J', |/'| = 2-^'-^, I J'l = 2-^'-^ + 2-^'-^,j' > k'. Applying the above discrete 
Calderon reproducing formula and the estimates in Corollary 3.5 yields for all (tt, v) e R, 

R'^I'xJ',j'>k' ' ' ' ' ' ' ' ' 





\r\ 








(1^1 


+ 


\u 






(1 J'l 


+ 1' 


V — 





R'=I'xJ',j'<k' 



A ^)^X 





\r\ 










+ 1 


\u 




xr\)(^+^) {\r\ 


+ 


\v 







where K, L are any positive integers which can be chosen by L, K > ^ — 1( for general n, m, K 
can be chosen by K > (n V m)(| — 1)), the constant C depends only on K, L and functions ip 
and 0, here xj' and j/j/ are any fixed points in J', respectively. 

Adding up all the terms with multiplying |/|| J| over RCQ, we obtain 
(5.1) J2 |/||^l'5i?<C J] J]|/'||J'|r(i?,i?')P(i?,i?')T^i?', 



RCQ, RCn R' 



where 



and 



-(^'^0 = (^a^)^-^(Ha^)^-i 



P(i?, R') 



if j' > k', and 



(1 _|_ distjl,!') y+K^i _|_ dist{J,J') y+x 

P(R R') — - 

^ ' ^ (^l I distjl,!') | dist{J,J') y+x 
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if / < k'. 

To estimate the right-hand side in the above inequahty (5.1), where we first consider R' = 
r X J', = 2-^'-^, \J'\ = 2-^'-^ + 2-^'-^,j' > k'. 

Define 

Jl»'^= IJ 3(2*/ X 2V) for z,£ > 0. 

Let be a collection of dyadic rectangles R' so that for z,£ > 1 

Bi^i = {R' = 1' X J', 3(2V X 2V) Pi n''^ and 3(2^-^/' x 2^" V) p| Q''^ = 0}, 

and 

Bo,£ = {R' = /' X J',3(/' X 2V)n^°'^ ^ and 3(/' x 2^-^00^°'^ = ^) ^ - 1' 

and 

Bi^o = {R' = 1' X J', 3(2V' X J') Pi 17*'° 7^ and 3(2*"^/' x J') f] 1]*'° = 0} for i > 1, 

and 

Bo,o = {R' :R' = l'x J', 3{I' X J') f] ^ 0}. 

We write 

To estimate the right-hand side of the above equality, we first consider the case when i — i = 0. 
Note that when R' e So,o, 3R' 0. For each integer /i > 1, let J^h = {R' = I' x J' e 
So,o, 1(3/' X 3J') nJ^°'°| > (2^)|3/' X 3J'|}. Let = Th\H-u and = [Jr^^v^ R' ■ Finally, 
assume that the right-hand side in (1.12) is finite, that is, for any open set f2 C i?^, 

\I\\J\TR<C\n\^-\ (5.2) 

R=IxJCQ. 

Since i?o,o = Uft,>i ^ft and for each R' e -Bo,0: P{R, R') < 1, thus, 

5^ Y,\i'\\J'HR,R')nR,R')TR' 

R'&Bo,o RQ^ 

^E E Ei^'ii^'w^o^ii' 

h>l R'CQ^ RCQ 

For each h > 1 and i?' C we decompose {i? : i? C f2} into 

Ao,o(i?') = {i? C n : dist(/,/') < |/| V |/'|, dist(J, J') < \ J\ V |J'|}; 
yli,^o(i?') = {RQfl: 2*'-^(|/| V |/'|) < dist(/,/') < 2*'(|/| V |/'|), dist(J, J') < \ J\ V |J'|}; 
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Ao,e'iR') = {RQfl: dist(/,/') < |/| V |/'|, 2^'-'^{\J\ V \J'\) < dist(J, J') < 2^'(|J| V |J'|)}; 
Ai,^t{R') = {RCQ: 2''-^{\I\ V |/'|) < dist(/,/') < 2*'(|/| V |/'|), 
2^'-\\J\ V \J'\) < dist(J, J') < 2^'{\J\ V I J'l)}, 

where i',£' > 1. 

Now we split ^ ^ \I'\\J'\riR, R')PiR, R')Tr> into 

h>l R'Cnn RCQ 

xr{R, R')P{R, R')Tr, =: + /2 + ^ + h. 
To estimate the term Ii, we only need to estimate r{R,R') since P{R,R') < 1 in 

ReAo,o{R') 

this case. 

Note that R e Ao,o(-R') implies 3Rf]3R' ^ 0. For such R, there are four cases: 

Case 1: |/'| > |/|, \ J'\ < |J|; Case 2: |/'| < |/|, \J'\ > |J|; Case 3: |/'| > |/|, \J'\ > |J|; Case 
4: |/'| < |/|, \J'\ < \J\. 

In each case, we can show Yr^Ao o '^^■^^ — C2~^^ by using a simple geometric argument 
similar to that of Chang-R. Fefferman [CF3]. This, together with (5.1), implies that Ii is bounded 

by 

J]2-'^^|nft|f-^ < c5]/ii-i2-'^(^-i+i)|nO'0|i-i < C\n\i-\ 

h>l h>l 

since \Qh\ < Ch2''\Q^^^\ and |O°'0| < C|0|. 

Thus it remains to estimate the term I^, since estimates of I2 and Is can be derived using the 
same techniques as in Ii and I4. This term is more complicated to estimate than term Ii. 

For each i' J' > 1, when R e A,/,^/ (J?'), we have P{R,R') < 2-*'(i+-^)2-^'(i+-^). Similar to 
estimating term /i, we only need to estimate the sum YreA , ^, '^(-^5 ^')- Note that R G A^/ £/(f?') 

implies that 3(2*'/ x 2^' J) n 3(2*'/' x 2^' J') ^ 0. We also split into four cases in estimating this 
sum. 

Case 1: |2*'/'| > |2»'/|, \2^' J'\ < \2^' J\. Then 



J^^^|3(2*'/' X 2^' J')| < |3(2*'/' X 2^' J') n 3(2*'/ x 2^' J)| 

< C2*'2^'|3i?'nl7°'°| <C2'2^'^^\3R'\ < C^^|3(2*'/' x 2^' J')|. 

Thus 1 2*' /'I = 2^*- 1+^1 2*'/ 1 for some n > 0. For each fixed n, the number of such 2* /'s must be 
< 2^ • 5. As for \2^'j\ = 2"*|2^' J'| for some m > 0, for each fixed m, 3 • 2^' Jn 3 • 2^' J' ^ implies 
that the number of such 2^ J' is less than 5. Thus 

J?£casel m,n>0 
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Similarly, we can handle 



Case 2: 


\2''l'\ 


< 


\2''l\, 


\2''J'\ 


> 


|20|, 


Case 3: 


\2''r\ 


> 


|2^7|, 


\2''J'\ 


> 


\2^'J\ and 


Case 4: 


\2^'r\ 


< 


|2^7|, 


\2''J'\ 


< 


\2''J\. 



Combining the four cases, we have YIr^A-, ^i{R') '"'i-^i ^C2 , which, together with the 
estimate of P(-R, i?'), implies that 

h<CY, E E 2-^^2-^'(i+^)2-^'(i+^)|/'||j'|T^,. 

h>l i',£'>l R'Qilh 

Hence /4 is bounded by 

h>l h>l 

since \Qh\ < Ch2''\n^^^\ and < C\Q\. 

Combining Ii, I2, I3 and /4, we have 

Now we consider 

E E E i^'ii-^'w^om^OT^i?'- 

i,e>i R'eBij flcn 
Note that for R' e 5^,^, 3(27' x 2^ J') n fi^'^ ^ 0. Let 

jrj^ = {i?' e Si,^ : |3(27' X 2^ J') n > ^|3(27' X 2^J')|}, 

and 



Since = IJ/i>i first estimate 



for some i,i,h> 1. 



Note that for each R' e V]l, 3{2W x 2^ J') n O^-^'^-i = 0. So for any C O, we have 
2^(|/| V |/'|) < dist(/, /') and 2^(| J| V | J'|) < dist(J, J'). We decompose {R : i? C 17} by 



Ai.,e{R') = {RGn: 2''-^ ■ 2\\I\ V |/'|) < dist(/, /') < 2'' ■ 2^(|/| V |/'|), 
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2^'-i • 2^(1 J| V I J'l) < dist( J, J') < 2^' • 2^(1 J| V I J'l)}, 
where i',£' > 1. Then we write 

J2 E l^'ll<^>(^'^')^'(i?,i?')T^i?'= EE E \I'\\J'\r{R,R')P{R,R')TR> 

Since P{R,R') < 2-^(i+^)2-^(i+^)2"*'(i+^)2-^'(i+-^) for R' e B,j and R G Ai,^e,{R'), r^ 
peating the same proof with Bq q replaced by Bi^i and a necessary modification yields 

J2 E \I'\\J'HR,R')P{R,R')Tr' < C2-^(i+^)2-^(i+^)2-^'(i+^)2-^'(i+^)x 
R'ev'/ ReAi,^^,{R') 

zl-i2^(i-i)£i-i2^(i-^)/.i-^2-'^(^-l+^)sup^U y \I'\\J'\Tr'. 



Adding over all i, £, i', /t > 1, we have 



^3tE E Ei^'ii'^>(^'^')W^')^^i?'<csup^l-^ 5^ iriij'iT«,. 

Similar estimates hold for 

E E E i^'ii-^x^'^on^'^o^fl' 

and 

E E E 

£>i R'eBo,i Rca 

which, after adding over all i,£ > 0, shows Theorem 1.21. We leave the details to the reader. 
Q.E.D. 

As a consequence of Theorem 1.14, it is easy to see that the space CMOp is well defined. 
Particularly, we have 

Corollary 5.1. We have 

^ R^supi EE E f{xi,yj)\^\I\\J\ 



CMO^ 



j k IxJCQ 



where I,J are dyadic cubes in R^,R^ with length 2 ^ ^,2 ^ -^-1-2 ^ ^, and xi,yj are any 
fixed points in I, J, respectively. 



Before we prove Theorem 1.16, we remark again that Theorem 1.16 with p = 1 in the one- 
parameter setting was proved in [FJ] on R"^ by use of the distribution inequalities. This method 
is difficult to apply to multi-parameter case. We will give a simpler and more constructive proof 
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which also gives a new proof of the result in [FJ]. Moreover, this constructive proof works also 
for other multi-parameter cases. 

Proof of Theorem 1.16: We first prove cP C (s^*)*. Applying the proof in Theorem 4.4, set 
and 

= {{x,y) eR^xR^: s{x,y) > T}. 

Denote 

Bi = {{IxJ) : |(/x J)nni| > i|/x J|,|(/x J)nl)i+i| < ^|/x J|}, 

where /, J are dyadic cubes in R^, R^, with side length 2~-'~^, and 2~^~^ + 2~^~^ , respectively. 
Suppose t = {tjxj} e and write 

I ^ SixjtixjW 
IxJ 

—\ sixjtixj\ 

i iIxJ)€Bi 

i (/xJ)GBi (IxJ)€Bi 

(5.2) <C||t||ep{$^|ar-^{ ^ k/xj|'}^}^ 

i {IxJ)el3i 

since if / x J e then 
and {t/xj} e yields 

{ixJ)eBi 

The same proof as in the claim of Theorem 4.4 implies 

J2 \sixj\^ <c2'''m 

{IxJ)€Bi 

Substituting the above term back to the last term in (5.2) gives C (s^)*. 

The proof of the converse is simple and is similar to one given in [FJ] for p = 1 in the 

one-parameter setting on i?". If ^ e (s^)*, then it is clear that £{s) = ^ sixjtjxj for some 

IxJ 

t — {tixj}- Now fix an open set O C i?" x R^ and let S be the sequence space of all s — {sjxj} 
such that I X J CQ. Finally, let be a measure on 5' so that the /i— measure of the "point" IxJ 
is — ^— -. Then, 

{ i^/xjp}^ = \\tixj\\p{s,df^) 

l^r /xjcn 
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, 1 - , 

= sup I j-j 2^ S/xJ*/xj| 

ll«IU2(S,dM)<l IxJCQ 

1 

< ||t||(sP)* sup \\SixJ iZtIIsp- 



IU^(S,dM) — 

By Holder's inequality, 

1 

P/X J 2 tWsP 



IxJCQ 
< 1, 



which shows HtHcp < [[^[[(s^)*- Q.E.D. 

In order to use Theorem 1.16 to show Theorem 1.17, we define a map 5' which takes / e 
(Sf)' to the sequence of coefficients {s/xj} = {\I\^\J\^'^j,k * fixi,yj)}, where I, J are cubes 
in R^,R^, with side length 2~^~^ ,2~^~^ + 2~^~^ ^ and xi,yj are any fixed points in /, J, 
respectively. For any sequence s = {s/xj}, we define the map T which takes s to T{s) = 

Y.Y.Y.Y.\^\^\J\^i^j,k{^^y)^ixJ^ where i^j^k are the same as in (1.9). 
j k J I 

The following result together with Theorem 1.16 will show theorem 1.17. 



Theorem 5.2:. The maps S : ^ and CMO^p c^, and T : ^ and CMO^p 
are hounded, and T o S is the identity on Hp and CMOp. 

Proof of Theorem 5.2: The boundedness of S on Hp and CMO^ follows directly from the 
Plancherel-Polya inequalities. Theorem 1.9 and Theorem 1.14. The boundedness of T follows 
from the same proofs in Theorem 1.9 and 1.14. To be precise, to see T is bounded from to 
Hp, let s = {s/xj}- Then, by Proposition 4.1, 

\\Tis)\\np < ^IKE E E E l^^-'^ * T{s){x, y)Wi{x)xj{y)}^\p. 

3 k J I 

By adapting a similar the proof in Theorem 1.9, we have for some < r < p 

\tpj,k * T{s){x,y)xi{x)xj{y)\^ 

=\Y.Y.Y.Y.\^'\\^'\^^^^*^rA-r){x,y)srxj^^^^^^^ 

j' k> J' r 

<c E 2H^-^'i^2H'=-'='i^{M,(EEi*^'x^ni^rvrVx.')'^}'(^,i/)x/(^)xj(y) 
k'<j' J' I' 

+ j2 2H^-^'i^2H'=-'='i^{M(EEi^^'x^Hi/rvr'xj-,o^}'(^,i/)x/(^)xj(i/)- 

k'>j' J' I' 
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Repeating the proof in Theorem 1.9 gives the boundedncss of T from to Hp. The similar proof 
given in the proof of Theorem 1.14 apphes to the boundedness of T from to CMO^. We leave 
the details to the reader. The discrete Calderon reproducing formula, Theorem 1.8, and Theorem 
1.14 show that T o is the identity on and CMO^. 

We are now ready to give the 

Proof of Theorem 1.17: If / e «Sf and g e CMO^, then the discrete Calderon reproducing 
formula. Theorem 1.16 and Theorem 5.2 imply 

l^gl ^\<f,9>\ = 

I ^ \I\\J\^R* f{xi,yj)i^R{g){xi,yj)\ <C\\f\\H^\\g\\cMO''J,■ 
R=IxJ 

Because Sp is dense in H^, this shows that the map ig =< f,g>, defined initially for f e Sp 
can be extended to a continuous linear functional on Hp and \\£g\\ < C\\g\\(jMQP^. 

Conversely, let £ G (Hp)* and set £i = £ oT, where T is defined as in Theorem 5.2. Then, 

by theorem 5.2, £i e (s^)*, so by Theorem 1.16, there exists t = {t/xj} such that £i{s) = 

Yl ^ixjtixJ for ah s = {sixj} and \\t\\cp ~ ||-^i|| < C\\£\\ because T is bounded. Again, by 
IxJ 

Theorem 1.16, £ = £ o T o S = £i o S. Hence, with f e Sp and g = J2 tixJ^Pnix — xj, y — yj), 

IxJ 

where, without loss the generality, we may assume ip is a, radial function, 

£(f)=£,(S(f)) =< Sif),t>=< f,g>, 

This proves £ = £g and, by Theorem 1.16, |b||cMO^ < C'll^llcf < C\\£g\\. Q.E.D. 

Proof of Theorem 1.18 As mentioned in section 4, since Hp is a subspace of L^, by the 
duality of Hp and BMOp, and the boundedness of flag singular integrals on Hp, one concludes 
that L°° is a subspace of BMOp, and flag singular integrals are bounded on BMOp and from 
L°° to BMOp. This shows Theorem 1.18. Q.E.D. 

6. Calderon-Zygmund decomposition and interpolation on flag Hardy spaces Hp{K^ x 
i?^): Proofs of Theorems 1.19 and 1.20 

The main purpose of this section is to derive a Calderon-Zygmund decomposition using func- 
tions in flag Hardy spaces. Furthermore, we will prove an interpolation theorem on H^{R^ x R^). 

We first recall that Chang and R. Fefferman established the following Calderon-Zygmund 
decomposition on the pure product domains i?^ x i?^ ([CF2]). 

Calderon-Zygmund Lemma: Let a > be given and / G L'p{R^), 1 < p < 2. Then 
we may write f = g + b where g e L^{R^) and b e H'^{R\ x R\) with \\g\\l < a^-^\\f\\l and 
Il^llifi(it2^xii2^) ^ C'q;-'^~p||/||^, where c is an absolute constant. 

We now prove the Calderon-Zygmund decomposition in the setting of flag Hardy spaces, 
namely we give the 
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Proof of Theorem 1.19 We first assume / e L2(i?^+"^) n x i?"*). Let a > and 

Qe = {{x,y) e i?^ X i?"^ : S{f){x,y) > a2^, where, as in Corollary 4.6, 



[ i,k I,J 

It has been shown in Corollary 4.6 that / e L^{R'^+'^) n H^iR"" x i?"^) then ^ \\S{f)\\p. 

In the following we take i? = / x J as all dyadic rectangles in i?" x R'^ with |/| = 2"-^"^, 
I J| = 2~^~^ + 2~^~^ ^ where j. A; are integers and N is large enough. 

Let 

7^o = |i? = / x J, such that \R n Oo| < ]^\R\ 

and for £ > 1 

= |i? = / X J, such that |i?nf)£_i| > but \Rr\VL(,\ < ^l^lj 

By the discrete Calderon reproducing formula in Theorem 4.4, 

f{x,y) = '^^\I\\J\^jk{x - xi,y -yj)(t)jk* {T^^if)) ixi,yj) 

j,k I, J 

= Y1 l^lWjkix - xi,y -yj)(f)jk* {T^\f)) {xi,yj) 

£>1 IxJeTZi 

+ Yl \I\\J\^jk{x-xi,y-yj)(f)jk*{T~'^{f)){xi,yj) 

IxJeTZo 

=b(x,y) + g{x,y) 
When pi > 1, using duality argument, it is easy to show 



\9\\p^<C\\\ Yl 

\(f>jk* {Tj^ '(/)) {xi,yj)\'^xixj 

M=IxJ€'Ro 



\Pi • 



Next, we estimate 11(71 when < < 1. Clearly, the duality argument will not work here. 
Nevertheless, we can estimate the Hp^ norm directly by using the discrete Calderon reproducing 
formula in Theorem 1.8. To this end, we note that 



\\9\\hI^ ^\\\YY1 K'^j'k' *g){xr,yj')\^xi'{x)xj'{y) \ Ikfi- 

Since 

('0j',fe' *9){xi',yj') = ^ |/||^| (iJj'k' *^jk^ {xr -xi,yj' -yj)(t>jk* (^i^^l/)) {xi,yj) 

ixJeTio 
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Repeating the same proof of Theorem 1.9, we have 



j',k'I',J' 



<c\\\ Yl \<l>jk*{TN\f))i^i,yj)Wixj 



ipi- 



This shows that for all < pi < oo 



I^IIh^i <c||<j Y \4>jk*{TN\f))i^i,yj)\'xixj 



2 

I Ipl- 



Claim 1: 



/ SP^{f){x,y)dxdy >C\\\ ^ \<Pjk * (^/, 2/j)I'x/XJ 

Js(f)(x,y)<a o_r., 



'S(/)(x,j/)<a U=/xJe7io 

This claim implies 



II^IU <C / SP'{f){x,y)dxdy 

JS{f){x,y)<a 

<CaP'-P I SP{f){x,y)dxdy 

Js(f)(x.v)<a 



'S{f){x,y)<a 

<CdF^ ^ll/lli?^(ii"xfl'")- 
To show Claim 1, we denote R = I x J & TZq. We choose < g < ^'i and note that 

/ SP^{f){x,y)dxdy 

Js{f){x,y)<a 



Sif){x,y)<c 



YY\^jk * (Ti,\f)) {xi,yj)fxi{x)xj{y) \ dxdy 

j,k I, J 



PI 

2 



dxdy 



PI 

{xuyj)\'^XRr\Q.i{x-,y)\ dxdy 



^jk * (^Ar^(/)) {xi,yj)\'^XRnQ^) {x^y))" 

Jr-xR^ [ Ue7^o 

PI 

>c/ I 5^ \<Pjk*{T^\f)){xi,yj)\\R{x,y)\ dxdy 



- ' 9 
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where in the last inequahty we have used the fact that Ifig fl (/ x J)| > ||/ x J| for / x J e TZq, 
and thus 

and in the second to the last inequality we have used the vector- valued Fefferman-Stein inequality 
for strong maximal functions 



1 1 

^ OO \r / OO \r 

Y.(Ms{fk)Y] \\p<c\\lj2\fkr] Hp 

Kk=l J \k=l ) 



with the exponents r = 2/q> 1 and p = pi/q > 1. Thus the claim follows. 
We now recall = {{x,y) e i?^ x i?"^ : Ms(xnJ > |}. 
Claim 2: For p2 < 1, 

II Yl \m\M^-^i,y-yj)<f>jk*{Ti:^\f))ixi,yj)\\^^.^,<c^^^ 

Claim 2 implies 

ig., <5^(2v^|fv;i 

l>l Js{f)ix,y)>a 

<Ca^'-^ ! S^{f){x,y)dxdy < Ca^^-^l |/| |^. . 

Js{f){x,y)>a 



To show Claim 2, again we have 

II mM^-xi,y-yMjk*{Ti^\f)){xi,yj)\\% 



P2 
F 



j'k' I', J' 



Y (^o'k' * (l>3kj {xr - xi,yj> - yj)(j)jk * {T^^{f)) {xi,yj] 



IxJeTZe 



2^ 2 



(/)) {xi,yj)\'^xiXj 

M=IxJeTZe 



IP2 



where we can use a similar argument in the proof of Theorem 1.19 to prove the last inequality. 
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oo 



However, 

>[ S{fnx,y)dxdy 

= L \Y.Y.\^^>'*iTN\f))i^i,yj)\'xi{x)xj{y)\ dxdy 

jQe-iX^e \^ j^k I,J J 



P2 

2 



P2 
2 



£2 
2 



P2 
2 



> 



/ J J] |(/.,fc*(T-H/))(x/,yj)px/(x)xj(y)i dxdy 



In the above string of inequalities, we have used the fact that for R gTZi we have 

1, , , , 1, , 

\Rnne_i\ > -\R\ and |i?nn^| < -\R\ 

and consequently R C fi^-i. Therefore |i?n > ^\R\. Thus the same argument applies 

here to conclude the last inequality above. Finally, since L^(i?"+"^) is dense in H^[R^ x i?"^). 
Theorem 6.1 is proved. Q.E.D. 

We are now ready to prove the interpolation theorem on Hardy spaces for all < p < oo. 

Proof of Theorem 1.20: Suppose that T is bounded from to L^^ and from to L^^. 
For any given A > and / e i?^, by the C alder on- Zygmund decomposition, 

f{x,y) ^ g{x,y) + b{x,y) 

with 

II^II^V, < CX^^-Wh^ and < CA^-^||/||^. . 

-Tip, 

Moreover, we have proved the estimates 

Ibll^V <C f S{fY-{x,y)dxdy 

" Js(f)(x.v)<a 



and 

|P2 



g., <C /" S{fY\x,y)dxdy 

^ Js{f){x,y)>a 
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y)| > ^}\da 



which imphes that 

poo 

\\Tf\\P =p / a^-'\ {{x,y) : \Tf{x,y)\ > A} \da 
Jo 

<pj^ aP-^\i^{x,y):\Tg{x,y)\>^yda+pj^ a^-'l i^{x,y) : \Tb{x,: 

POO P POO P 

<p / a^~^ / S{fY^{x,y)dxdyda + p I a^~^ I S{fY'^{x,y)dxdyda 

Jo JsU){x,y)<a Jo J S{f){x,y)>oc 

<c\\f\\V, 

Thus, 

I|T/||p<C||/||h^ 

for any p2 < p < pi. Hence, T is bounded from to L^. 

To prove the second assertion that T is bounded on Hp for p2 < p < Pi, for any given A > 
and / e Hp, by the Calderon-Zygmund decomposition again 

\{{x,y):\g{Tf){x,y)\>a}\ 
<\ {{x,y) : \9{Tg){x,y)\ > |} | + | : b(T6)(a;,y)| > |} | 

<Ca-P' I \Tg\\P^r^ + Ca-P' \ \Tb\ 
<Ca-P'\\g\fj^,,+Ca-P'\\b\f^, 



tP2 



<Ca-P^ I S{f)P^{x,y)dxdy + Ca-P^ I S{fY^{x,y)dxdy 

Js{f){x,y)<a Js{f){x,y)>a 

which, as above, shows that ||T/||jjp < C||gr(TF)||p < C||/||7jp for any p2 < p < pi. Q.E.D. 
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